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Abstract 

In this work we modify the wave-corpuscle mechanics for elementary charges intro- 
duced by us recently. This modification is designed to better describe electromagnetic 
(EM) phenomena at atomic scales. It includes a modification of the concept of the 
classical EM field and a new model for the elementary charge which we call a balanced 
charge (b-charge). A b-charge does not interact with itself electromagnetically, and 
every b-charge possesses its own elementary EM field. The EM energy is naturally 
partitioned as the interaction energy of pairs of different b-charges. We construct EM 
theory of b-charges (BEM) based on a relativistic Lagrangian with the following prop- 
erties: (i) b-charges interact only through their elementary EM potentials and fields; 
(ii) the field equations for the elementary EM fields are exactly the Maxwell equations 
with proper currents; (iii) a free charge moves uniformly preserving up to the Lorentz 
contraction its shape; (iv) the Newton equations with the Lorentz forces hold approxi- 
mately when charges are well separated and move with non-relativistic velocities. The 
BEM theory can be characterized as neoclassical one which covers the macroscopic as 
well as the atomic spatial scales, it describes EM phenomena at atomic scale differently 
than the classical EM theory. It yields in macroscopic regimes the Newton equations 
with Lorentz forces for centers of well separated charges moving with nonrelativistic 
velocities. Applied to atomic scales it yields a hydrogen atom model with a frequency 
spectrum matching the same for the Schrodinger model with any desired accuracy. 
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1 Introduction 

It is well recognized that the classical electromagnetic (CEM) theory formulated in the form of 
Maxwell-Lorentz equations provides an excellent description of electromagnetic phenomena 
at the macroscopic length scales. It also well known that CEM theory is inadequate in 
explaining electromagnetic (EM) phenomena at the atomic scales including spectroscopic 
data of the hydrogen atom (HA). 

We develop here an EM theory which accounts for all classical EM phenomena at the 
macroscopic scales as well at least some EM phenomena at the atomic scale including the 
HA spectral lines. Our theory is classical, though we can apply it to some phenomena (HA 
energy spectrum in particular) at spatial scales compared with Bohr radius, that is spatial 
scales of order 0.1 nm, and our theory produces the same type of results as the quantum 
mechanics which is commonly used at such scales. We think that expanding the classical 
theory down to smallest possible spatial scales is important because the classical description 
of physical systems allows (at least in principle) more details compared with the probabilistic 
quantum-mechanical description. Since our neoclassical description applied to HA does not 
generate contradictions at atomic scales of order 0.1 nm (as CEM theory did), we expect it to 
be applicable to describe details of electromagnetic processes at nanometer scales. Note that 
at such scales the electron cannot be considered a point or a charged ball and we have to use 
a complete description of an electron which is presented in our model as a wave- corpuscle. 

When attempting to change the CEM theory we want: (i) to stay on solid ground of the 
Lagrangian mechanics and the relativity principle; (ii) to recover in this new EM theory all 
well established experimental facts described by the CEM theory. The foundational pillars 
of the CEM theory - the Maxwell equations and the Lorentz force expression - remain to 
be key elements in proposed here new EM theory, and before we proceed with the new 
developments let us briefly recall the CEM theory fundamentals. First of all, the EM fields 
driven by prescribed currents in vacuum are described by the Maxwell equations 



where E and B are respectively electric field and the magnetic induction, and p (t, x), J (t, x) 
are respectively prescribed charge and current densities. In particular, when the fields are 
generated by point charges, the sources p, J are written in the form 



-^ + Vx E = 0, V ■ B = 0, 



c dt 
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where q l is charge value of the £-th point charge, r £ and is v £ (t) = are respectively its 
position and the velocity, and 6 is the Dirac delta-function. Second, for a given EM field the 
motion of every point charge in the field is determined from the equation 



d 



dt [mV {t)] = <f E (t, r< (*)) + V (t) x B (t, / (t)) 



(4) 



where m is the £-th point charge mass and the right hand side of (j3J) is the Lorentz force. The 
CEM theory essentially treats three types of problems: (i) studies of EM fields for prescribed 
charge and current densities described by the Maxwell equations (PQ), (EJ); (ii) the motion of 
charges in a prescribed external field; (iii) interaction of charges and their EM fields. The 
classical Maxwell-Lorentz system though very successful in describing many EM phenomena 
has well known problems including infinite self-energy which are discussed widely in the liter- 
ature. As to important for the CEM theory concept of point charges or particles F. Rohrlich 
writes, |Rohl Section 2.1, p. 9; Section 6.1, p. 123]: "... the classical theory of charged particles 
as first conceived by Lorentz emerged as a hybrid theory of particles and fields: charged par- 
ticles are interacting via an electromagnetic field", "Macroscopic Maxwell electrodynamics 
knows only charge distributions. In this theory electrostatic charge is derived from a charge 
density (linear, surface, or volume density). The concept of charge as an aggregate of el- 
ementary charged particles is foreign to it." There is also a fundamental thermodynamical 
problem for the classical Maxwell theory at atomic scale namely the lack of an "elementary 
process of absorbtion" as it was formulated by A. Einsteins in his seminal paper |Einl909b] . 
He wrote there: " The fundamental property of the oscillation theory that engenders these 
difficulties seems to me the following. In the kinetic theory of molecules, for every process 
in which only a few elementary particles participate (e.g., molecular collisions), the inverse 
process also exists. But that is not the case for the elementary processes of radiation. ... 
The elementary process of emission is not invertible. In this, I believe, our oscillation theory 
does not hit the mark. Newton's emission theory of light seems to contain more truth with 
respect to this point than the oscillation theory since, first of all, the energy given to a light 
particle is not scattered over infinite space, but remains available for an elementary process 
of absorption." In particular, the proposed here theory does have an elementary process 
of asbroption as A. Einstein suggested. We call this process "negative radiation" and its 
essense is that the EM energy is moving toward two closely located elementary charges if 
they oscillate with the same frequency but oposite phases. 

The proposed here EM is based on a new concept for elementary charge which we call 
balanced charge or b-charge for short, and we refer the theory itself as balanced electromag- 
netic (BEM) theory. A key element of the BEM theory is a concept of an elementary EM 
field assigned to every single b-charge. A single b-charge is described by a pair [ijj,A^), 
where ip is its wave function and = (<£>, A) is its 4-vector elementary potential with the 
corresponding elementary EM field defined by the familiar formula = d^A u — d v N 1 . So, 
b-charge is a field over 4- dimensional space-time continuum for which the wave function ip 
represents its matter properties and the elementary potential A^ mediates its EM interactions 
with all other b-charges. Importantly, (i) all internal forces of a b-charge are exclusively of 
non- electromagnetic origin; (ii) every b-charge is a source of its elementary EM field which 
represents force exerted by this charge on any other b-charge but not upon itself. The later 
allows to view a single b-charge as truly elementary one with respect to the electromagnetic 
interactions. 

An idea to introduce an extended charge instead of the point one is not new, and the 
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most known models for it are the Abraham rigid charge model and the Lorentz relativistically 
covariant model. These models are studied and advanced in many papers, see |Jac| Sections 
16], |Kie2] . |Pealj . |Roht Sections 2, 6], |Schwin] . |Spohn| , |Yaghjian] . In contrast to those 
models, here and in BabFigl] , |BabFig2] we do not prescribe to an elementary charge a cer- 



tain geometry, but instead the elementary charge has a wave function governed by a nonlinear 
Klein-Gordon or a nonlinear Schrodinger equation in the relativistic and nonrelativistic cases 
respectively. An idea to eliminate self-interaction is also, of course, not new. The latest to 
our best knowledge attempt to have this feature in the electrodynamics is due J. Wheeler 
and R. Feynman, |WFlj . | WF2] . but the EM theory proposed here is very different from it. 

The BEM theory is constructed based on a relativistic Lagrangian with the following 
properties: (i) b-charges interact only through their elementary EM potentials and fields; 
(ii) the field equations for the elementary EM fields are exactly the Maxwell equations with 
proper conserved currents; (iii) a free charge moves uniformly preserving up to the Lorentz 
contraction its shape; (iv) the Newton equations with the Lorentz forces hold approximately 
when charges are well separated and move with non-relativistic velocities. Since an over- 
whelming number of EM phenomena are explained within the CEM theory by the Maxwell 
equations and the Lorentz forces the BEM theory is equally successful in explaining the same 
phenomena. 

A system of N elementary charges in the BEM theory is modeled by N pairs (-0 , A^) , 1 < 
I < N, and every b-charge (ip e , A e>M ^ is naturally assigned via the Lagrangian its elementary 
conserved 4-current J eu . The later property provides additional justification for calling the 
elementary field F ^ v electromagnetic. The classical (total) EM field = (<£>, A) is recovered 
in this theory as the sum of all elementary EM fields, namely 

<p= E a= Al > ( 5 ) 

i<e<N i<e<N 

but, importantly, this total field is not an independent entity with its degrees of freedom. 
Notice then that since in the BEM theory there is no EM self-interaction the action on £-th 
charge by EM fields of other charges is described by a field A^ e = (<^^, A^J which is the 
total field (<p,A) "balanced" by the removal from it the self- interaction, namely 

<p& = XV > A +l = J2 (6) 

Use of EM fields similar to ones in is, of course, not a discovery and they can be found in 
many textbooks, but the BEM theory goes further than that and removes the elementary EM 
self-actions for the CEM Lagrangian and consequently the elementary self-energies from the 
classical EM energy-momentum tensor. Only with this removal of the elementary self-actions 
from the classical EM Lagrangian one gets the field equations in the form of the elementary 
Maxwell equations d^F efJjU = ^J iv for the elementary EM F e ^ u fields with the elementary 
conserved 4-currents J lv . 

In Section [2] we study differences and similarities between the CEM and BEM theories. 
The differences and similarities stem from the differences and similarities between their La- 
grangians and with consequent differences and similarities in the energy-momentum tensor 
definition. Since the both theories are based on the Maxwell equations they have exactly 
the same EM fields for prescribed currents but their total actions and the energy-momentum 
tensors differ with consequent differences in radiation phenomena. We show below that the 
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CEM theory is a limit of the BEM theory. That can be already seen by comparing formulas 
([5]) and for a large N since these large sums differ just by one term. The classical charge 
can be treated within the BEM theory as a cluster of many b-charges and it is shown, in 
particular, that such a classical charge does interact with itself electromagnetically and its 
interaction with other classical charges can be described approximately through the inter- 
action with the single effective EM field. Notice also that the classical single EM field is 
defined by ([5] as the sum of all elementary EM fields, but since the energy is a quadratic 
function of fields, the total single classical EM field has energy equal to the combined energy 
of interaction between all the pairs of elementary EM fields only approximately. Most of EM 
phenomena at the macroscopic scale can be described based on the single EM field with 
astounding precision and the relative difference with the BEM theory since the inverse of the 
Avogadro constant, that is 10 -23 , is very small. But the differences between the CEM and 
BEM theories become more pronounced for smaller systems with fewer b-charges. The BEM 
theory predictions can significantly deviate from those of the CEM theory in the following 
situations: (i) there are just a few b-charges which are in close proximity; (ii) there is a 
large but highly coherent system of b-charges similar to those collective, coherent systems 
(superconducting ring, laser and more) described by C. Mead |Mead| p. 5] in his "collective 
electrodynamics". In the BEM theory the EM energy is the energy of EM interaction of 
pairs of b-charges and consequently it is naturally partioned between the pairs of charges. 
We show that for every pair of b-charges their EM interaction energy satisfies elementary 
energy-momentum conservations governed by the relevant Lorentz force densities. 

In Section[3]we study particle-like properties of b-charges (t/> £ , A***) which are captured by 
the concept of wave-corpuscle similarly to the wave- corpuscle mechanics (WCM), |BabFigl|, 



BabFig2 . A key ingredient providing for particle-like behavior of a b-charge are self- 
interaction nonlinearities G a ^ where size parameter a = a e determines the size of a free 
particle which though small is non-zero. The shape of the wave-corpuscle tp is intimately 
related to the the nonlinearity G a / via the rest charge equation - a nonlinear Klein-Gordon 
equation- which in non-relativistic case turns into a nonlinear Schrodinger equation. We also 
derive in this section the Newton equations with the Lorentz forces as an approximation 
when charges are well separated and move with non-relativistic velocities. 

In Section we provide a detailed sketch of the hydrogen atom model. At atomic scales 
when there are just a few elementary charges in a close proximity the BEM theory differs 
significantly from the CEM theory. It yields, in particular, the BEM hydrogen atom model 
with a frequency spectrum matching the same for the Schrodinger hydrogen atom with any 
desired accuracy. The difference between the two HA models depends on on the size of the 
free electron as a parameter in the self- interaction nonlinearity G a , and an analysis suggests 
the size of a free electron to be of about 100 times of the Bohr radius. 



2 Relativistic theory 

In this section we introduce BEM Lagrangian for a system of bounced charges and derive 
the corresponding field equations, currents, energy-momentum tensors (EnMT) for system 
components and more. An important objective of this section is to study the EM interactions 
in the BEM theory and find out common and different features of the BEM theory and the 
CEM theory. 
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2.1 Lagrangian, field equations and currents 

Let us consider a system of elementary charges A^) , 1 < £ < N. In what follows the 
i-th charge potential A 1 * and its EM field F ttaJ = A lv — A 1 * completely account for its 
action upon all other charges £' 7^ I. Consequently, the action upon £-th charge by all other 
charges is described by £-th exterior potential A^ and its EM field F? v defined by 

A^ = J2 Ai '^ 4 M =(^> A y> F l f = Y,F^ v . (7) 

We also introduce for the total potential and the corresponding total EM field by 
the following formulas 

1<£<N 1<£<N 



We furnish now the system of iV b-charges with the following Lagrangian 

£({^,v&},{^,v&M*) = E £M^,^S) + £bem, (9) 



KKN 



J^Vjr p^v pi 
^BEM = ^CEM — £e, ^CEM = — , £ c = _ / J — , (10) 



167T ^— ' 167T 

\<1<N 

where L l is the Lagrangian of the £-th bare charge, and the covariant derivatives are defined 
by the following formulas 

^ = aV, ^; = § e ^ e *, (ii) 

\<f A 1 * ia e A e>M 

where ip* is complex conjugate to ip 

Observe that EM part £bem of the Lagrangian C according to f fTUj) is obtained by the 
removal from the classical EM Lagrangian £cem all self-interaction contributions C c of the 
elementary EM fields and it can be recast as 

rpifMU rpi' -pi^v -pi 

= - E ^-E^- m 

{£,£'}:£'^£ 1<£<N 



The "bare" charge Lagrangians L are defined by exactly same expressions as in BabFigl 
BabFig2|, namely 

L e ^ ^S) = £j ~ «W - G e } , (13) 

where (i) G is a nonlinear self-interaction function of the £-th charge described below; (ii) 
m > is the charge mass; (iii) q e is the value of the charge; (iv) x > is a constant similar 
to the Planck constant h = and 

f f f 1 

f ur m c , m c , . 

/c = — = , ur = . (14) 

c X X 
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The system Lagrangian £ defined by (IT2"l) - (IT4"j) is manifestly Lorentz and gauge invariant with 
respect to the gauge transformations of the first kind (11181) . The gauge invariance via the 



Noether's theorem allows to introduce elementary conserved currents, |BabFigi] , |BabFig2 

.q e ( dL 



J 



-1 — 



x \M V 



dU_ 



4 



*£ 



dU 



dA 



with the conservation law 

dvJ* = 0, V + V • 3 l = 0, J» = (p £ c, J £ ) . 



(15) 



(16) 



The Euler-Lagrange field equations for the above Lagrangian C are (i) elementary wave 
equations 



+ K t* + Q 1 ' f Ut\ 2 ) = o, = d> + 



together with the conjugate equation for ip 



= ^ 



(17) 



;is) 



and (ii) the Maxwell equations for the elementary EM fields 



c 



with the familiar vector form 



V • E* 



VxE' + -d t W = 0, 

c 



V • W = 0, 



Att . 



(19) 

(20) 
(21) 



VxB 1 - 

c c 

We will refer to the equations ( I19p . (|2"0"|) as the elementary Maxwell equations and to the field 
equations ( TT7]) - (|T9|) as field equations for b-charges. Using ( fill]) -( TT5"]) we obtain the following 
representation for the £-th elementary current J £u = (cp £ , J e> ) 



J 



fx 



fx \v 



Im 



2m e 
ip £ XC 



(22) 



or its vector form 



e \ it- 
q \ip 



P 9 

]£ = f \ip | 



x Im —r + 9 

~4? ~ 



Xlm 



(23) 
(24) 
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The Maxwell equations for elementary EM fields (fl9|) combined with the equalities (JHJ)-© 
readily imply that total and exterior fields also satisfy the Maxwell equations 

(25) 
(26) 



d^ u = 


c 


where J v = 


1<£<N 




An 

TV 

c ^ 


where = 





Observe that: (i) the field equations ([19]) for the elementary EM fields are exactly the 
Maxwell equations with the corresponding elementary currents; (ii) every elementary wave 
equation ffTTj) - ([T9]) indicates that the £-th charge is driven by its exterior potential A*£ in- 
dicating that there is no self-interaction. We can see also from equation (|25p that the total 
EM field satisfies the Maxwell equations for the total currents as in the CEM theory. 

It is instructive to see how the BEM system Lagrangian £ defined by ([9l)- (TP2l) can be 
obtained by a modification of a similar WCM Lagrangian involving the classical single EM 
field introduced in , BabFigl , BabFig2|. The modification consists of two actions: (i) 



alteration of the covariant derivatives fill I) which removes the self-action by using exterior 
potentials A 1 ^ instead of the same total field potential A 11 ; (ii) subtracting from the classical 
action of the total EM field the sum C e of the classical actions of the individual EM fields 
as in (J9j). It is the two described actions combined yield an EM theory with the Maxwell 
equations and the Lorentz force densities as its exact components. The natural alteration 
of the covariant derivatives alone is not sufficient since then the field equations will be not 
exactly the Maxwell equations nor there will be force densities described exactly by the 
Lorentz formula. Importantly, an additional action - the subtraction of the Lagrangian 
component C e defined in (Q - is necessary to have the Maxwell equations and the Lorentz 
forces as fundamentally exact parts of the new EM theory. One can see already significant 
differences brought into the EM theory by many elementary EM fields. It is evident from 
the formula (I12p for £bem that the EM action is the sum of elementary EM actions £g EM 
associated with all pairs {£, £'} : £' ^ i of elementary charges and every elementary action 
£g EM depends on the fields F l ^ u and F e ^ v only. Continue this line we observe that though 
the EM Lagrangian £bem according to formula ffl2l is a simple and natural summatory 
function of the elementary EM fields, it can not be reduced exactly to any function of a 
single EM field as in the CEM Lagrangian. 

In what follows we often use a vector form of the system Lagrangian £ defined by 0, 
ffTOj). that is 

^(Wii.{(^A / )}y=W (27) 




S 



where 



C 



CEM 



£bem 

1 

8^ 



{21 



- E 



KKN 



d t A 



d f A l 



(V x A) 2 
- (V x A e ) 



The corresponding field Euler-Lagrange field equations are, first of all, the Maxwell equations 
for all elementary EM potentials ip e , A e 



-d t A l + V<p e = -4 np e , 



Vx(VxA') + -9 t ( -d t A £ + V<A = —J* 
v ' c \c / c 



l,...,iV, 



(29) 
(30) 



where the charge densities and currents are defined by (l23~i) . ( J24l) . and, second of all, the 



equations for the wave functions tp 1 in the form of the nonlinear Klein- Gordon equations 



0. 



l,...,iV, 



(31) 



and similar equations for the complex conjugate variables if)*. Note that equations ( 13TT) for 
ip are coupled with the equations for EM potentials via the covariant derivatives ((9]). If we 
choose for all elementary potentials the Lorentz gauge 



-<V + V ■ A E = 0, 

c 

then equations (1291) - (1301 turn into the wave equations 

vy - \dW = -W 

\d^A £ - V 2 A e = — 3 £ , 
c z c 

as in the CEM theory. 



l,...,iV, 



(32) 

(33) 
(34) 



Importantly, as in the WCM, BabFigl , BabFig2 , the nonlinearities G are determined 



based on the single £-th charge equations (j29l - (jHTl) under the assumptions that it is isolated 
and is at rest. Namely, let us consider a single b-charge, set iV = 1 in (|29|) - (|3T|) and simplify 
notations ip = ip , ip = A = A £ , G e = G. Then according to ([9]) the covariant derivatives 

are df = dt, V = V and the Lagrangian for a single b-charge is 



2m 



d t ip 



G 



(35) 



Observe that the Lagrangian Lq above does not depend on the potentials <p, A since there is 
no self-interaction. Though we can still find the potentials based on the elementary Maxwell 
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equations (|33|) . they have no role to play and carry no energy. Let us consider now the 



rest state of the b-charge which we set as in the WCM, |BabFigl , BabFig2 , to be of the 
form 

o 

^ (t, x) = e" la;o V (x) , co = = ck , (36) 

lp (t, x) = Cp (x) , A (t, x) = 0, 

where ip (|x|) and lp = lp (|x|) are real-valued radial functions. Substituting the ip, lo and A 
defined by the relations (136]) into the field equations (l3Tj) (l3Tj) we obtain the following rest 
charge equations: 



-V V + G' 



-V 2 <£ = 4vr 



■0 



^ = 0, (37) 

(38) 



The quantities ip and lp are fundamental for our theory and we refer to them, respectively, as 
form factor and form factor potential. The equation (I3TI) signifies a complete balance of the 
two forces acting upon the resting charge: (i) internal elastic deformation force associated 
with the term —Aip; (ii) internal nonlinear self-interaction of the charge associated with 

the term G' (^^>\j V*. We refer to the equation (|37|) . which establishes an explicit relation 

between the form factor ip and the self-interaction nonlinearity G, as the charge equilibrium 
equation. Hence, if the form factor %p is given we can find from the equilibrium equation 
(1371) the self-interaction nonlinearity G which exactly produces this factor under assumptions 
that ip (r) is a nonnegative, monotonically decaying and sufficiently smooth function of r > 0. 

Now we pick the form factor %p considering it as the model parameter and then the nonlinear 
self interaction function G is determined based on the charge equilibrium equation (37|]. As 
in the WCM, IBabFigl] , \BabFig 2], such a determination of the nonlinearity is a key feature 
of our approach: it allows to choose the form factor ip and then to determine matching 
self-interaction nonlinearity G rather than to deal with solving a nontrivial nonlinear partial 
differential equation. 

To explicitly integrate the size of the resting b-charge into its model we introduce size 
parameter a > into G = G a through its derivative as follows 

G' a (s) = a- 2 G[ (a 3 s) , where G' (s) = d s G (s) . (39) 

In the following Section [3~2l we give examples of the nonlinearity G and discuss its properties. 



2.2 Energy-Momentum tensors and Lorentz forces 

The conservations laws, particularly the energy-momentum conservation, play an important 
role in the physics of EM phenomena. A general source of conservation laws in Lagrangian 
theories is the Noether's theorem, [Goldl Section 13.7], which yields canonically conservations 
laws based on the Lagrangian symmetries, that is its invariance with respect to continuous 
groups of transformations. The conservation laws are not independent equations, and they 
hold only if the fields satisfy the field Euler-Lagrange equations. The energy-momentum and 
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charge conservations are two important laws in any EM theory and they have special signif- 
icance in our studies for several reasons. First of all, the conservation of energy-momentum 
describes their transport in the space and is directly related to the point charge approxima- 
tions. Second of all, the Lorentz force density expression which the one of the most important 
components of any EM theory arises in the energy-momentum conservation equations 



and not in the field Maxwell equations. Third of all, the elementary charge and energy- 
momentum conservation equations f|T6|) . fl58|) account for charges individuality. In fact, the 
BEM Lagrangian is invariant with respect to a wider group of elementary gauge transforma- 
tions (11201) . (11211) . and it is due to this symmetry the "two-way" representation f|T5|) holds 
for every elementary current: one that involves the differentiation with respect to the wave 
functions derivatives and another one that involves the differentiation with the respect to 
the elementary potentials. It is due to this two-way representation the conserved Noether's 
elementary current is exactly the source current in the corresponding elementary Maxwell 
equations ffT9|) . f l20|) . 

To get further insight into b-charges properties we need to find the symmetric EnMT 
of the system Lagrangian C defined by (Tr2"j) - (1T4"|) making use of a general method used in 



BabFigl , BabFig2 . The method yields the following representation 



'j'fiu ^ ^ rpifiu _|_ ^[iv 'zzliv ^ ^££' flf (40) 

1<^<7V I'^-l 

where the individual EnMT T ifll/ of the bare £-th charge and the EnMT components 1 E ie '^ u 
for the EM fields are as follows 

t*» = |^ + - rtf = (4i) 



X 



2 



47T \ 75 4 7€ 



?' ± £. (42) 



Notice that the expression ( )4Ti) is the same as in (BabFigl | , |BabFig2"] with the only difference 
that the covariant derivatives here are defined by equalities (1111) . To emphasize notationally 
the new meaning of the EnMT of the EM fields we use for it the symbol S MV and we continue 
to use the symbol G MI/ for the classical EnMT of the the EM field as in formula (I212p . As 
to the expression (1421) for the individual components of the symmetric EnMT it is obtained 
from the representation and evidently is similar to the expression (I212p for the classical 
EnMT . An alternative and more elementary way to derive the expression f|42|) is based 
on the expression (12101) for the canonical EnMT as it is done in |Jac[ Section 12.10]. 
Indeed, in our case the canonical EnMT takes the form 

E u ^ = i- + <T^f- , £' + L (43) 

47T l07T 

The identity d v A 1 '^ = -F^ + d^A l ' u allows to recast the expression (T431) into 

= — ( g^FtF^" + -g^F'F e 'A - —g^F^A^. 
4n \ K 4 K J 47r fS 



(44) 
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Then using the Maxwell equations ( TT9|) we can recast the second term as 

_ Lgmpi^A^ = # (_±. g w F ^ A ^ _ (45) 

The term —J^A 1 ' 11 combined with the canonical EnMT T liiV yields the the symmetric 
and hence 

1<£<AT 

l<i<N I'^l J 

The later by the standard argument implies the conservation law 

d^V lu = (47) 

for the symmetric EnMT T^ v defined by the expressions ( 140]) - (142]) . 

Observe that using representations ( I42p we can also recast expression (T4T)]) for EnMT 
as follows 

e'^e i<e<N 

where 



" 47T 



" 47T 



(49) 



The relation between the classical and the new EnMTs and H Miy in view of their repre- 
sentations (HDD, (S2J) and ( 12T2D is as follows 



I'jJL 



gez'fw = J p^ v J _ ^ p^^ v ^i = (5Q) 

\1<£<JV / l<^<iV 

1<£<N 

Evidently the formula fl50|) reads that the new expression of the EnMT equals to the classical 
EnMT of the total EM field minus the sum of the the the classical EnMTs for the elementary 
EM fields. 

Based on the formula (|42p for the components S for £' ^ i of the symmetric EnMT we 
can represent their entries in terms of the fields and B^, namely 

8n 

■'' " " " 1- ' 
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1 

3^ 



E^Ef + BfBf - -S^ (E e ■ E e + B e ■ B e ' 



W Cj 
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(52) 



(53) 



which are evidently similar to the same entries (I213l) - (l215p for the classical EM theory. In 
particular, we have the following expression for the Poynting vectors similar to (12181) 



= c 2 g\ 



cE £ x W 

47T 



(54) 



To find expressions for the EM force densities acting upon b-charges let us examine their 
conservation laws. Combining the Maxwell field equations ffl~9]) . ( 126]) . an elementary identity 



g^F^ - frF^ = d v F^, for F>* = d"A^ - &A>", 
and the antisymmetry of the EM field tensors F ttu/ we obtain for £' ^ £: 

1 ji F e '^ + — --F 1 d v F 1 '^ + -ff (F e F e '^\ 



(55) 



(56) 



Let us introduce an elementary EM field interaction energy for the pair {£,£'} of b-charges 
which in view of (156]) satisfies the following elementary energy-momentum conservation laws 



(5* 



with the right-hand side being the negative of the sum of the corresponding Lorentz force 
density. 

It readily follows from relations ( 157]) and (149]) that the interaction EnMT has the 

following representation in terms of the classical EnMT 



^ {£,£'} Liv _ qiiv ^plfiv _|_ pl'uvj _ qui/ ^pt^v^ _ q[iv ^pi'ni/j 



(59) 



The representation (159]) is evidently similar to (150]) . 

Notice that in view of (|5T]) . ( 152]) and (154]) the entries of the tensor «w l/^ are as follows 



w 



{£,?} 



;{V}oo 



• W + B* 



c 9i 



{£,£'} _ ^{i,e'}oi _ ^{i,e'}io 



4tt 

E e x B e> + E t> x B * 



(60) 



An 



g{£,£'} 



E £ x B e ' +E f 'xB 

47T 
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The expressions ( 160]) and (16T1) can be alternatively derived from their classical counterparts 
(l2T3j) . d2HD and (12181) based on relation ([59} 

The relations ( 140]) together with (1581) readily imply the total energy-momentum conser- 
vation law 



(62) 



l<€<iV 



with the sum of the negative of the Lorentz forces in the right-hand side. 

Using expression (14 lj) for and the same transformations as in the similar case con- 



sidered in BabFig2 we obtain the following elementary conservation laws 



1 



f TP* 



(63) 



Observe that the expression on the right-hand side of the above equality is the Lorentz force 
density 4- vector acting upon the 4- vector current Jj by its exterior field F, , and the same 
vectors with the minus sign arise in the conservations laws for the EM fields (1581) and (1621) . 
Notice also that the natural partition (1401) - (T4"2"l) of the system EnMT into charges EnMTs 
T^ v and the EnMTs H^'" for pairs {£, £'} of interacting elementary EM fields produces 
exactly Lorentz force densities in the elementary conservation laws (|58|) . ( 162|) and ( )63|) . This 
provides a solid justification for the energy-momentum partition pUjl - pS]) . and it seems there 
is no simple way to alter it when preserving the Lorentz forces expression. 

The vector form of the elementary conservation law ( l5"8j) is similar to the CEM theory 
conservation law 12201 |Jacl| Section 6.8], namely 



d t w^ + V • S<^> 



(64) 



dg 



dt 



where ^ is the Lorentz force density satisfying 



(65) 



ft + f!\ ft = P l Ef + -(J f xB { 



(66) 



Remark 1 The conservation law (64\ ) can be easily derived from the classical energy conser- 
vation law (the Poynting theorem) as follows. Let us consider the classical energy conservation 
law associated with the 3 Maxwell equations / TJPj) : (i) fori-th EM fields, (ii) for i'-th EM field 
and (Hi) for their sum. Then subtracting from the conservation law for the sum of the fields 
the sum of conservation laws for l-th and £'-the fields we obtain exactly the conservation law 



2.3 Elementary EM fields for prescribed elementary currents 

To study the properties of the elementary EM fields for balanced changes and the energy- 
momentum transfer in the space-time it is instructive to consider a situation of prescribed 
currents for b-charges similarly to the same for the case the CEM theory Lagrangian ( I207j) . 
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In the BEM theory the corresponding Lagrangian £emj is based on the representations (J9]) 
and is of the form 



C 



BEMJ 



167T 167T C 

l<£<Ar KKN 



Based on the representation (fl2|) it can be alternatively written as 

£ BEMJ = - E - \ £ 4-4? (68) 

{£,£'}:£' y^e KKN 
pliiv pi 1 

4^ 16tt c ^ ^ 

1<£<7V l<£<iV 

The field equations for the Lagrangian £emj can be obtained directly from the variational 
principle. Indeed, taking into account the antisymmetry of every field tensor F £liu , the fact 
every pair £, £' with £' ^ £ appear two times in the Lagrangian representation fl68l) and the 
usual assumption of the decay of all the fields at infinity we obtain 



5 [c BEM jdx= j (^F^-hA5A%dx 

"* KKN ^ ^ 



(69) 



Since the variables 5A^ U , 1 < £ < N can vary independently (see relations (I122p and (I123p ) 
the requirement for the above variation to vanish yield the field equations which are exactly 
the same Maxwell equations for the elementary EM fields as in (fl9|) . (|20|) . ( 120]) . 

Importantly, the system of elementary EM fields with prescribed currents satisfy the 
elementary conservation laws (|6%|) . (|65|) which can also be directly derived from the Maxwell 
equations (12~TT) . Indeed using the general vector identity (\272\\ and the Maxwell equations for 
the indices £ and £' we obtain 

V • (E< x B*') = B £ ' ■ (V x E e ) -E f -(vx B*') = (70) 

= -B e ' ■ -d t B e -E e -( -d t E £ '+— J"' 
c \c c 

Adding then to the identity (ITU]) a similar one obtained from it by swapping indices £ and £' 
we obtain the elementary energy conservation (j6"%|) . Similar direct derivation is possible for 
the elementary momentum conservation ( 165]) . 

Observe also that according to the BEM conservation laws ( T55|) the energy and the momen- 
tum are assigned not to the elementary EM fields by themselves but rather to their interacting 
pairs. This is a noticeable change compared to the CEM theory where the EM field has an 
energy and a momentum on its own. 



2.3.1 Dipole elementary currents 

For simplicity sake let us assume all prescribed elementary currents to be in the form of ideal 
electric dipoles. Recall that an ideal electric dipole source concentrated at a point xq, with 
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the charge and current densities Jd and p d , is defined as follows, [vBla| Section 7.10, (7.151), 
Appendix 8], 

J (t, x) = d t [p (t) 5 (x - x )] =p(t)8(x- x ) , (71) 
p (t, x) = -V • [p (t) 5 (x - x )] = -p (t) ■ V6 (x - x ) 

where 5 (x — xo) is the Dirac delta-function, and p (t) and dtp (t) satisfy 

Pi (t) = J x' jP {t^) dx', j = 1,2,3, p = {pi,p 2 ,Pz), (72) 

p(t) = d t p(t) = y J(t,x')dx'. (73) 
It readily follows from (I7T1) (1237P that the potentials of the ideal electric dipole are 

p(*o) 



p(t,x) = 0, A(t,x) 



c|R| 



where 



R = x — xq, R 



R 

|RJ 



tn = t 



x - Xq 



(74) 



(75) 



Then applying the Jefimenko and the Panofsky-Phillips formulas (1241 j) and (12421) for the ideal 
electric dipole sources p (t, x) and J (t, x) defined by formulas (1711) we obtain the following 
formulas for the EM field 



E(t,x) 



3 R-p(to) R-p(to 



IRI 



(76) 



3 (p (t ) • RJ R - P (to) (p(t )xR)xR 



B(t,x) 



c|Rr 

'p(to) p(*o) 



c 2 IRI 



x R, where p = d^p. 



(77) 



cR 2 c 2 R 

When deriving formula (1751) we used a vector identity (12691) . 

In a simpler case when the dipole function p (t) is time harmonic and complex valued of 
the form 

P (*) = Pu,e- iart , p (t ) = P.e fc l R l~^, where k = -, (78) 

c 

the general formulas ( 176|) . ( 1771) yield the well known formulas for the ideal electric dipole 
fields, [Jacll Section 9.2], 



£,2 e fc|R|-iu;t 

E(t,x) = — x 



(79) 



x < p x R x R 



3 R- p R- p 



FIR!" k R 



B(t,x) 



R 



fc R 



50) 
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implying in radiation zone k |R| ^> 1 the following asymptotic formulas 



E (t, x) = fcV |R| ~ iw< 



R x x R 



IRI 



1 + 



A; IRI 



^1) 



B(t,x) = fcV |R| - iw * 



R x p a 
IRI 



1 + 



fcR 



For the case of a multi-frequency dipole function p (t) we can introduce 

min { — \ > 0. 



weAp I c 



(82) 



(83) 



Then radiation components decaying as |R| 1 dominate for /c min |R| 3> 1 in formulas (1761) (T77|) 
implying the following asymptotic expressions for the radiation fields, see also [Gril Section 
11.1.4], 



E(t,x) 



p'(t ) x R) x R 



c 2 IRI 



1 + 



B(t,x) 



p' (to) x R 
c 2 |R| 



1 + 



h ■ R I 
""min | xl '| 

for A; m i n IRI > 1. 



Energy flux for a system of arbitrary elementary dipoles Using expression (18"4"|) for 

the EM fields and the vector identities (12701) and (12711) we consequently obtain the following 
formulas for the energy flux S u and the corresponding total powers P u and P^ e ' £ ^ radiated 
through a sphere centered at the dipole location 



4+ 



-W x 



(85) 



p e (t ) ■ r (to) - (V (t ) ■ r) (r • r (to; 



4ttc 3 |Rr 
p^(to)xR)-(y (to)xR 



R 



4vrc 3 IRI 



^R, 



P" = [ S* da = ±p< (t ) • p* (t ) , (86) 

J|x| = |R| 3C 3 
oC 

Let us assume now that the dipole functions p £ (t) depend on time t almost periodically as 
in Section 16.41 Then the representation (181)]) together with relations ( 16"Uj) . ( 175]) and (I265p 
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consequently imply the following formulas for the time-averaged radiated powers 

<^-pO-((V-r) (ft-p* 



4ttc 3 IRI 2 



R = (87) 



"((> xR,) 


• 






)] 


47TC 3 


R 


2 



R 



5_ £ c 4 Re{(pf ; xR)(p^xR)} 



8ttc 3 |R 



we A ^nA ./ 
p p 



(88) 



weA ^nA 
p p* 



4(p'.p<") „ 



toeA i(~)A pl 
p 

where A p t and A pl i are respectively the frequency spectra of p e (t) and p e (t). It readily 
follows from the formula then that i/ie time- averaged radiated power (p{W}\ can take 
any real value: negative, zero or positive. Indeed, according to the formulas 



[P^j vanishes if the frequency spectra A p t and A pl i don't have any common frequencies, 



TP* 

i.e. 



pm) = 4 <p'W c f (')) =Oitvnv = . (90) 

In particular, the relation $9U\) shows that if the both l-th and £'-th b-charges are monochro- 
matic of different frequencies then the time-averaged radiated power is exactly zero or in other 
words there is no radiation. The formulas f l89|) . f!265j) readily imply 



p^) = tmym = _i e w , k , » > if ^ (t) = p < (t) , (91) 

wga i 
p 

= _i»WW> = - A E - 4 | P £,| 2 < o if „«' (*) = -p* (t ) . (92) 



weA £ 
p 



Evidently the relation KM) describes a situation when for a given pair {£, £'} the time-averaged 
radiated power (p{W}^ is negative, that is the radiated energy propagates with the speed of 
light toward the source rather than away from it. 

Radiated power of a system of identical dipoles A basis for simple comparison of 
radiative properties in the BEM and CEM theories let us consider a system of N > 1 b- 
charges described by one and the same dipole moment p (t). Then in view of f[9Tj) we have 
for the system the following time-average radiated power 

<iw = J2( p{ '- n )= < 93 > 

weAp 
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In the CEM theory a similar system of N identical dipoles p has an effective dipole moment 
Np. Substituting this number for p in the expression (1551) we get 

2(p* (t)-p' (g) 1 

weAp 

Relating representations (193]) and (1941) we readily obtain 

(Pbem) = (l - (Pcem) • (95) 

2.3.2 Elementary currents for point charges 

It is curious to see how the BEM theory describes fields and radiations phenomena when 
the elementary currents are caused by a system of point charges as, for instance, in the 
Rutherford atom model. In the case of point charges we use the acceleration fields E a (t, x) 
and B a (£, x) defined by formulas (I253p . (I254p to find the time-averaged radiation power. 
Consequently, for every pair {£, £'} of b-charges we have 

s te = x g^' and s et = _£.eJ x B{. (96) 

47T ' 47T 

Suppose that one of the two b-charges, say £-th charge, is at rest or moves uniformly implying 

that is f3 =0. Then it follows from formulas (12531) . (I254p that E^ = and B e a = implying 
that S u and S vanish, namely 



if j3 = or fa =0. (97) 



Suppose that the position functions r e (t) and r e ' (t) of the corresponding b-charges are almost 
periodic functions as described in Section 16.41 and that mod (V and mod (r^ (t)) have 
no common frequencies. Then in this case the time-averaged flux ^S^') is exactly zero, i.e. 



if mod (/ (t)) p| mod (r e (t)\ = 0. (K 



Now let us compare this estimate with a similar computation in the framework of the CEM 
theory with the total field E = Ef + E e a , B = Bf + B^ or in the framework of the BEM 
theory with two prescribed charges and a test charge. In the framework of the BEM theory 
if we consider action of the system of two charges onto a test charge or onto a distant system 
of test charges which does not noticebly affect the radiation of two charges, the field of two 
charges turns into an external field given by the same formula E = Ef + E£, B = B* + B£. 
The radiation power at large distances generated by two charges in accordance with ( I252p is 
expressed by the formula 

S = ± (Ef + Ef) x (Bi + Bf ) . (99) 

If £-th charge is at rest B e a = and E^ = and we obtain the leading part and its time 
average 

S = |-EfxBf, (S) = ± (Ef x B f) , (100) 
which differs from (19"8"|) and ( 19"7|) . 
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Remark 2 Note that formulas / TPS)) and (9l\ ) describe the radiation of EM interaction energy 
between the two charges, and it is not the same as the EM radiation of the fields which act 
on test charges. One may ask why we are interested in the interaction energy? The answer 
is that this is the energy which is derived from the Lagrangian and the interaction energy is 
conserved when the system is closed. Of course the system where the motion of sources is 
prescribed is not closed, and there must be an external source of energy to provide for the 
prescribed currents and for the radiation which acts onto test charges. 

Let us consider in conclusion the case when the both charges move slowly and hence 
|/3 £ | , 1/3 I <C 1. Notice that that relations (1255 j) defining the point charge EM fields are 
similar to the dipole expressions (j84"|) and can be obtained from them by a substitution 
p £ (t) — > cqir (t). Consequently, expressions for the energy flux and the radiated power 
for the pair of slowly moving charges can be readily obtained from the dipole expressions 
(!85|) - (l85|) by a substitution p £ (t) — > cqv e (t) with an additional correcting asymptotic factor 
5 e = (1 + (\/3 £ \)). In particular, 



?2 ¥ (t) ■ V (t) - (V (t) ■ R.) (R • ¥' (t) 



47TC 



IRI 



K5 



a' 



q 2 


"(V (t) x r) 


■ (V x R) 




Aire 47rc 3 


|R 


2 



R<r , where 
5 U ' = {l + 0{\(3*\))(l + 0(\[3 e '\)). 

2.4 Comparison with the classical EM theory 

We would like to show here that the CEM theory can be viewed as a limit case of the 
BEM theory Particularly, we are interested in seeing: (i) how the single classical EM field is 
modeled by elementary EM fields in BEM theory; (ii) how classical EM phenomena including 
the radiation are modeled within the BEM theory. 

As an example of relation between the classical EM field and elementary EM fields we 
consider clusters of many tightly bound identical b-charges, every cluster is labeled by index £. 
Namely, for every I we introduce Ne b-charges having identical wave functions and elementary 
EM potentials with charges = Npq 1 : 



(t,s) 



(t,8) 



A (£,s)n = A (i,s)n where I < £ < N, l<s<N e 



(101) 



and identical EM fields F l » v = F^ s ^ u and currents J e > u = J^' s > . Now we compare these 
fields with classical EM fields defined by ([8]) with the same currents by comparing the 
electromagnetic parts of the Lagrangians in BEM and CEM theories. We have for classical 
theory 



C 



CEM 



167T ^ 



1 



16tt 



16tt 

V NpNpF^FL - 



V NeNpF^Ff' 



16tt 



N'F^Fl. 
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For the BEM theory (jl~2l) takes the form 



C 



BEM 



167T ^{(e,s),(i',s'):(i',s')^(e,s)} 

— V N e ,N £ F e ^F e ' u - - — V N e (N e -l) F e " u F e u 



The difference between £cem and £bem can be attributed to to interactions inside every 
cluster. In particular, we have for the both theories 

Ccemi = NfF^F^, £ B em£ = N e (N t - 1) F^F^ 
readily implying the following expression for the relative difference 

£bem^/£cem<? — 1 = (102) 

The difference (11021) evidently becomes small as number Ng of particles in the cluster becomes 
large. 

To summarize, the CEM theory with its single EM field can be formally derived from the 
BEM theory as a limit obtained by binding together N# identical particles with Ng — > oo. 

2.4.1 Lagrangians for clusters of charges 

Here we make a comparison of two Lagrangians for N clusters of charges with Ng particles 
within every cluster which are formally derived from WCM and BEM theories respectively. 
To this end we consider first a Lagrangian as in WCM theory with Ni + ... + Nn charges and 
fixed nonlinearities G 1 ^ 

To this end we consider first a Lagrangian as in WCM theory with N\ + .. . + Nn charges 
and a fixed nonlinearity 

c ({^}, {t$*>}) = E E L " (v^,^) - (los) 

l<i<N l<s e <N e 

and to introduce clusters we impose additional constraints as in (11011) . namely 

lP = l£, q l ' s = q^ s = l,...,N e . (104) 

A cluster in the equilibrium is described by the real valued functions ip and tp e,s which 
satisfy the following system of equations: 

-A^ = 4.£Wl-^Wf , (105) 
-A/- + (2 - + G' w (|**|>) i* = 0. (106) 



Note that the WCM theory developed in BabFigl , BabFig2 coincides with the case where 



every cluster contains exactly one charge. In this case the charge equilibrium equation has 
the form 

-A^ = W (l-^) tf\ (107) 
-A^^( 2 -^' + G ;(|^)^0. (108) 



21 



Comparing with f 1 1 5 j) and (I106P and setting 

<fi = N e tf, qi = N e q e , ml = N?m e (109) 

^' 5 = ^ = V4, s = l,...,N e , 
we find then that (p^ satisfies the following equation 

-A^ = 4rr4(l-#iU*, (HO) 

which has exactly the same form as (11071) in the WCM theory. The charge normalization 
condition is also fulfilled. Equations ( 11051) and (11081) can be rewritten in the form 

- a + 'Ml 6 - MV + <?. fw'i 2 ) *' = 0. tin) 



This equation has the form of (1 1 6 f) for = ip w if we set 

Xw = Nfa. (112) 

Observe that relations (11141) - (11091) readily imply the following identities for the model con- 
stants 

mi,c m c 



(113) 

If we introduce Lagrangian 

({VO> {VO) with (i) constants defined by (dUD, ([112) and 

(ii) EM potential defined by 

^ = 2 ^ = 2 A w , ^ = 1 < £ < JV, (114) 



then the Euler-Lagrange Equations for Lagrangian C ( {^' Se) }, {^ l ' H) } ) defined by (fT03|) 



with restrictions f ll04j) are equivalent to the Euler-Lagrange Equations for C w . Hence, we 
can conclude that introduction of clusters of charges is equivalent to a proper rescaling of the 
constants of the WCM theory. Note that the quadratic dependence of m m on the number Ni 
of charges in the cluster is natural in the relativistic theory since the energy of interactions 
depends on the number of particles quadratically, and relativistic mass of the cluster is 
proportional to its energy. 

Now let us relate the above treatment with the BEM theory. To make the comparison 
more transparent we consider a generalized version of the BEM theory, namely the theory 
where the basic equilibrium object is a cluster of N# b-charges, Ng > 1. The BEM theory for 
elementary b-charges corresponds to the case Nf = l for all I in this generalized setting. The 

Lagrangian £bemc {\^^ St) }-, {V^'^H f° r clusters is given by an obvious modification of (jUJ) 

and the variables are subjected to additional constraints as in (I104p . The equilibrium state 

for a generalized cluster of n = Ni identical b-charges is denoted by f bn ,^' n , s = 1, ...,n. 
The equilibrium condition takes the following form similar to f l 1 5 [) and ( I106P 

N e / Ls ° I 

1 9bnVb,« 



s=2 \ m b,„ c2 



A VV + r^^b n 2 TJ-; V'b.n + G b,n [ Wb,n\ ) A,n = °" ( 116 ) 

Xb,n V m b\n C J V ' 7 
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Note that the only difference between fl 1 5 [) and fl 1 6 [) is that the first term with s = 1 in 

the sum is omitted. Now ipb2i'Pb2 is the basis for comparison and once again we fix the 
nonlinearity G' hn = G' b2 — G' h and set similarly to (11091) . (11121) : 

<fi,n = ( n ~ 1 ) V&.2, <lln = (n-l) q{ 2 , (H7) 
m i,n = in- if mj )2 , Xb,n = (n - I) 2 Xb,2- 

Let us introduce Lagrangian £ B em,h ({V'b.nHVb.n}) with so defined constants. Observe that 
the Euler-Lagrange equations for £bemc ({V^'^}; {^f^^} ) with restrictions fll04j) are equiv- 
alent to the Euler-Lagrange equations for jCbem,™ ({V^nHv^n})- Obviously if all > 2 then 
£bem,h. has the same form as £ w . In addition to that, if 

G b,2 = G w, 9b,2 = 1 , m b,2 = m , Xb,2 = X, 

then the relative difference of coefficients of Lagrangians £bem,« and C w with the same 
n > 2 is of order -. For example, (a^, — ql n ) /q^ = l/n. Notice that the case 77 = 1 (which is 
the primary one considered in this article) is special, and there is a non-vanishing difference 

between ^ b2 ,^ b2 and the fundamental equilibrium V'biJ^bi = ^ > *P determined by ( 1371) . 

m- 

Thus, based on the above analysis we conclude that the WCM theory can be considered 
as an approximation for the generalized BEM theory if the WCM charge is identified with a 
cluster of a large number of b-charges. 

2.5 New EM features of the BEM theory 

In this section we discuss new EM features of BEM theory not presented in the CEM theory. 
Looking at the b-charge Lagrangian £ with its field equations and conserved currents in 
Section 12.11 we can already see new important features of the theory of interacting charges 
and EM fields. The first striking feature of the new theory compare to the classical one is 
that the single EM field as independent entity is no more, instead we have elementary EM 
potentials A 1 ^ and fields F l ^ v defined by the classical formulas and satisfying individually 
Maxwell equations / TTPj) . Though we still can naturally define the total potential A 11 and the 
corresponding total EM field J- 7 ^ as the sum of individual potentials and fields by formulas 
©-(HD, and though the total EM field J=^ satisfies the Maxwell equations (J25D this total 
field is not an independent entity. 

With all that said about the new status of the total EM field its physical significance 
as the field sensed by a small test charge remains valid in the BEM theory. Indeed let 
us look at how different charges sense each other via the field equations (JT7j)-(fT9l). It is 
evident from the equations and the formulas ( 122]) for individual currents that the entire 
evolution of £-th charge is completely determined by a single quantity, its exterior potential 
A e ^, and, consequently, charges (ip e , sense each other exclusively by their potentials 
A e>M . In other words, the state of every b-charge (ip , A e>1 ) and its current J 1 ^ are completely 
determined by the action upon it of all other charges potentials via the exterior potential 
A e £ and every charge A efJl ) acts upon other charges. So in this theory as we can see it is 
the individual A 1 ^ potentials as components of respective b-charges (■?//, A 1 ^ that facilitate 
the EM interaction between charges. This is in contrast to the classical theory as well to 
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our own WCM theory in BabFigl , BabFig2 in which it is a single and independent EM 



field interacting with every charge facilitates all EM interactions between different charges 
■0 . In particular, in the BEM theory by its very set up there is no any EM self-interaction as 
in the classical theory and our theory in |BabFigl , BabFig2|, where such interaction exists 
because there is only a single EM field. 

Based on the analysis in Section |2~41 one may expect noticeable differences between CEM 
theory and BEM theory for example when 1/Ng in (I102p is not small. These differences can 
become more pronounced when C e in fl28|) is comparable with the classical EM Lagrangian 
£cem not to mention when if it is a dominant term in £bem- 

An important signature of the BEM theory differentiating it from the CEM theory is a 
mechanism of negative radiation for certain prescribed currents, i.e. a situation when the EM 
energy propagates with the speed of light toward the current source rather than away from it 
as we have shown in Section [231 This mechanism can conceivably work for a limited time in a 
system of several bound charges, such as an atom or a molecule, resulting in effective energy 
gain coming from matching energy loss of b-charges outside of this system. Such energy 
transfer is completely accounted for by interacting elementary EM fields as the components 
of the involved b-charges including external ones which might be represented effectively by 
an external EM field. The very possibility of the negative radiation indicates a significant 
difference between the EM energy transport at the elementary atomic scale from the same 
at the macroscopic scale. For the later most of the time one would observe well known 
classical charge behavior including the EM self-interaction and exclusively normal (positive) 
radiation of the EM energy away from the source. In other words the BEM theory allows for 
a differentiation between macroscopic and atomic scales at the level of EM fields alone, and 
that makes the concept of a b-charge with its elementary EM field to be truly elementary. 
This is in noticeable contrast with some of the classical charge models which are based on the 
concept of a single EM field, and we would like to quote here F. Rohrich, |Roh| Section 6-1]: 
"This was the problem faced by Abraham, Lorentz, and Poincare in the first few years of 
this century. The most obvious model of a charged particle is a sphere carrying a spherically 
symmetrical charge wave function. While such a model is meant (and was indeed proposed) 
as a picture of a charged elementary particle (an electron, for example), it is obvious that it 
is basically a macroscopic charged body, only much smaller. There is nothing "elementary" 
about it." 

Notice also that as we can see from the BEM expression (|60|) for the interaction energy 
density for a pair of b-charges the interaction energy density can be negative. The later is 
analogous to the negative sign of the electrostatic energy for two classical point charges when 
one of them is positive and another one is negative, and the interaction energy is defined to 
be zero when charges are separated by infinite distance. This way to calibrate the interaction 
energy in the CEM theory is in line with defining the interaction energy as work done to 
assemble a system of charges from the state when they don't interact, that is when they are 
separated by infinite distance. 

Let us look at the gauge properties of the BEM system Lagrangian C defined by Q, ffTUl) . 
Recall that gauge transformation of the first or the second kind (known also as respectively 
global and local gauge transformation) are described respectively by the following formulas, 
[PauRFl (17), (23a), (23b)], [Wen! Section 11, (11.4)] 

ip £ — > e 17 V £ , — > e" 17 V £ * ; where 7 £ is any real constant, (H8) 

„ iq^ A(x) „ „ iq^ X(x) „ 

ip £ ->• e — ~ip e , ip e * ->• e-i^-V , Atl -> A " + (119) 
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The significance of individual EM fields in the new theory is manifested in the existence of 
a new much large group of gauge transformations than described by (j!18p - (jll9p . We have 
already pointed out that the system Lagrangian C is invariant with respect to the gauge 
transformation of the first kind (11181) and used that to construct an instrumental for the 
theory conserved individual currents J lfi defined in (TTUT) . We introduce now a new gauge 
transformation which we call elementary gauge transformations gauge or transformations of 
the third kind: 

„ ig* A (x) „ „ lq" X (x) „ 

->• e —tfj £ , <tf* -> e^V*, (121) 

where functions X e (x) , 1 < £ < N, are independent real- valued scalar functions of x and 

A (x) = X " ( x ) > ^ ( x ) = X ( x ) ~ X " ( x ) > ( 122 ) 

1<£<N 

It is readily follows from (j!22p that 

\ e = \ - = - \\ where A = ^ X' = (N - 1) A, (123) 

i<e<N 

implying the independence of functions A (x) , 1 < £ < N. A straightforward examination 
shows that for iV > 2 the system Lagrangian L defined by (fl2"]) - (fll|) is invariant with respect 
to the gauge transformation of the third kind. As we have already pointed out it is due to the 
gauge invariance with respect to elementary gauge transformations (I120p . (11211) the two-way 
representation (fT5|) holds for the elementary currents that accounts for an important fact that 
the conserved Noether's elementary current is exactly the source current in the corresponding 
elementary Maxwell equations ( TT9l) . (1201) . 

Another feature of the BEM theory is the new expressions (HHl) . (H9l) and (1501) of the total 
EM energy and the entire EM field EnMT compare with the classical expressions (I212I) - (I214I) . 
Let us take a look at the representation fl50l) which relates the EM EnMT to the classical 
one. This identity shows that the EnMT can be effectively obtained by the removal from 
the classical EnMT for the total EM field the sum of all classical EnMT's for individual EM 
fields. So if our intention was to remove the EM self-interaction in a consistent way keeping 
the Lagrangian structure it is perfectly accomplished by the BEM theory. 



3 Non-relativistic dynamics of localized charges 

In this section based on our relativistic model we introduce a non-relativistic model for 
the case where charges move slowly compared with speed of light. First we introduce field 
equations and the Lagrangian and briefly describe their derivation from the relativistic model. 

Using frequency- shifting substitution (I3"6T) with more general ip = ip u (t, x) which depends 
on (t, x) we observe that the second time derivative in the Klein-Gordon equation (13TT) can 
be written in the form 

- = h (ft + 7^) 2 - a y + ^) € + 46 (is*) 
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where 

*£ = Xy. (125) 

To get a non-relativistic approximation, we neglect in (j!24p the term with the factor -7 
and substitute -2i^- (d t + ^<p^\ ifa + K&t instead of the term -^dfdfipl in As a 

result, the nonlinear Klein-Gordon equation (13 lj) is approximated by the following nonlinear 
Schrodinger equation 

^ + ^ <y) 2 ^ - ^ ^ - ^ = °' (i26) 

with V given by ( fill) , where for notational simplicity we write in place of iff u . Since 
magnetic fields generated by moving charges also have coefficient - and are small for small 
velocities, they can be also neglected, and consequently we preserve only the external mag- 
netic fields and replace V by the covariant gradient 

VL = V - (127) 

Notice that in the nonrelativistic case the £-th b-charge is described by a pair if)*, ip e where 
the elementary EM field is represented only by the scalar electric potential (p . The non- 
relativistic model has the the Lagrangian 

3vr 



£ (>> WL > WL) = l -^r + E 11 W ^ *) > where ( 128 ) 
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£ [tf*d*f - VW] - £7 { I V„^| + G l (W) 

Iv^l 2 



and 



¥> = £y. (129) 

Importantly, to total potential ^ defined by fl 129[) is not an independent entity but just the 
sum of the elementary potentials <p l , and if> e * a variable complex conjugate to tp e . 

The Euler-Lagrange equations for charge densities ifr take the form similar to (11261) . 
namely 

ix9rt4 ' = (^*) 2/ + qt y* + ^ ^ + ^ ki' d^O ^ <i3o) 

where ip^ e is given by f l 1 2 5 1) and every potential can be determined from the equation 

vV £ = -W|/| 2 , e = i,...,N. (131) 

The solution of the above equation f)13ip is given by the Green's formula 

V (t,x) = q / r dy. (132) 

Jr3 \y — x i 
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The properties and examples of the nonlinearities G £ a are provided in the following 13.21 As 
the result of charge conservation the norms HV^H remain constant on solutions of (I130p and 
choose those constants to 1. In other words we impose the following charge normalization 
condition: 

U e \\ 2 =f U*| 2 dx = l, t>0, £=1,...,N. (133) 

A motivation for this particular normalization is based on the formula f ll32[) and requirement 
that the elementary potential (p is asymptotically the Coulomb's potential q / |x| for large 
|x|. 

Recall that in the classical electrodynamics the evolution of a point charge q of a mass m 
and and position vector r (t) in an external electromagnetic (EM) field in the non-relativistic 
case is governed by the Newton's equation @. In our model a charge is described by a wave 
function ip (t, x) with dynamics described by a system of nonlinear Schrodinger equations 
(NLS) fll30p coupled through corresponding electric potentials. Nevertheless we show below 
that in macroscopic regimes we can derive from the field equations as an approximation the 
Newton's equations for centers T£ (t) of localized waves defined by 

r g (t) = / x|^ £ (t,x)| 2 (ix. 

These center adequately describe the positions of wave functions ifj e (i, x) when they are 
localized. To introduce localized solutions we use the nonlinearity which depends on a size 
parameter a > which, in turn, determines the spatial scale of a the charge when it is at the 
rest state. We consider below the macroscopic dynamics, with a macroscopic spatial scale 
R 3> a. We prove in Section H] that in the case of the non-relativistic field equations when 
a — > the centers of the interacting charges converge to solutions of the Newton's equations 
with the Lorentz forces if ip e remain localized. We also provide examples of exact solutions of 
the field equations in the form of accelerating solitons for which the localization assumption 
holds. 



3.1 Single charge 

In the case of a single charge the Lagrangian and the field equations are obtained by setting 
N = 1 in (I128p . (I130p . (I13ip and this case evidently <£>^ = 0. In particular, for a single charge 
without the external field we have (p 1 ^ = 0, (p ex = 0, A cx = and we obtain the following 
equations for the rest state ip with dtip = 0: 

? VV + + K]'(|V| 2 ) = 0, (134) 



2m* r ^ exr 2m e 

VV = -4:irq\ilj\ 2 . 

Obviously, the equations coincide with ( 1381) . and therefore the rest solutions of relativistic 
and non-relativistic equations coincide, as it should be expected in the case of zero velocity. 

Let us consider now a single charge, omitting the index £, in an external EM field and set 
N = 1 in f )130p . (113 ip . For a special class of external fields we present explicit solutions to the 
field equations (I130p . fl 1 3 1 [) for a single charge in the form of wave- corpuscles (accelerating 
solitons). We assume here for simplicity a purely electric external EM field, i.e. when A ex = 0, 
E ex (t, x) = — V<^ cx (t, x) (see |Ba bFigl] , |BabFig2 for a similar exact solution with non-zero 
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magnetic field). For the purely electric external field the field equations (I130p . (113 ip take the 
form 

= + qy e J> + ^G' a (\ip\ 2 ) V, (135) 

VV = -4vrg|^| 2 . (136) 
We define then wave- corpuscle ip, ip by the following formula: 

ip (t, x) = e [S/x ^ } S = mv • (x - r) + s p (t) , (137) 
$ = V>(|x-r|), p = £(|x-r|), r = r(i). 

In the above formula ip is the form factor satisfying ( j!4ip . <p is a radial function determined 
from f)13ip . and we refer to the function r (t) as wave-corpuscle center. Since ip is center- 
symmetric, this definition agrees with more general definition (I163p . 

Suppose that ip cx ((t, x)) is a continuous function which is linear with respect to x. Then 
ip defined by ( j 13TI) provide an exact solution to ( j 135j) . provided that r (t) is determined from 
the equation 

d 2 r (t) 

m^^ = qE ex (t,r), (138) 
and v (t) , s p (t) are determined by formulas 

v = ^, s p = jf f^f- - q Vcx (t, r (f ))) dt>. (139) 

The verification of the fact that (I137p , (I138p and f )139p determine an exact solution is straight- 



forward and details can be found in [BabFigl|, BabFig2 or BabFig3 



Note that in a simpler case when the external fields (p ex and A cx vanish, a simpler solution 
of fll35p -f fT36l) is provided by (I137P with r (t) = r + vt with constant velocity v. In this case 
the wave- corpuscle solution (I137P of the field equations fll35p -f fT36|) can be obtained from the 
rest solution ip,(p by certain Galilean-gauge transformations. Solutions of a similar form are 
known in the theory of nonlinear Schrodinger equations, see |Sul] and references therein. For 
the particular case of the logarithmic nonlinearity solutions of the form (I137P were found in 
[BiaJ in the form of accelerating gaussons. The exponential factor in wave-corpuscle solution 
of the form (1 1 3 T[) can be identified with the de Broglie wave, for details see |BabFigl 



BabFig2 . 



Remark 3 The construction of the solution does not depend on a particular form of 

the nonlinearity G' = G' a as long as §3l\) is satisfied. It is uniform with respect to a > 0, and 
the dependence on a in is only through ip (|x — r|) = a^ 3 ^ 2 ip 1 (a -1 |x — r|). Obviously, 

if ip (t, x) is defined by \13T^ then \ip (t, x)| 2 — > 5 (x — r) as a — > 0. 

Remark 4 In \BabFigiy , IBabFigj^j we introduced a non-relativistic model where equations 
for a steady states coincide with l[134\ )- We derived there nonrelativistic equations from the 
corresponding relativistic once not only assuming that velocities are vanishingly small, but 
also assuming that the Sommerf eld's fine structure constant is vanishingly small. If we would 
not neglect there the terms with Sommerf eld' s fine structure constant the steady states would 
satisfy more involved equations, see \BabFigI^ for details. We also would like to note that 
though ip would be defined differently, equations U38\) would be the same. In the present paper 
the rest state equations ( |i^| ) are simpler even if the Sommerf eld's fine structure constant is 
not replaced by zero. 
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3.2 Determination of Nonlinearity, 



As we have already mentioned, the nonlinear self interaction function G is determined from 
the charge equilibrium equation (137)) based on the form factor (ground state) -0. Important 
features of our nonlinearity include: (i) the boundedness or slow subcritical growth of its 
derivative G' (s) for s — > oo with consequent boundedness from below of the wave energy; 
(ii) slightly singular behavior about s = 0, that is for small wave amplitudes. 

In this section we consider construction of the function G, study its properties and provide 
examples for which the construction of G is carried out explicitly. Throughout this section 
we have 

■0, -0 > and hence |0| = ■0. 

We introduce explicitly the dependence of the free ground state 0> on the size parameter 
a > through the following representation of the function ip (r) 

(r) = i, a (r) = a" 3 ^ (a"V) , (140) 

where ip l (r) is a twice continuously differentiate function of the dimensionless parameter 
r > 0, and, as a consequence of f)133j) . the function ■0 O (r) satisfies the charge normalization 
condition for every a > 0. Obviously, definition (I14(jp is consistent with (157)) and ( 155)) . The 
size parameter a naturally has the dimension of length. A properly defined spatial size of ip a , 
based, for instance, on the variance, is proportional to a with a coefficient depending on 
The charge equilibrium equation (157)) can be written in the following form: 

V 2 ^ = G' a (#) Vv (141) 

The function ip a (r) is assumed to be a smooth positive monotonically decreasing function 
of r > which is square integrable with weight r 2 and we assume it to satisfy the charge 
normalization condition of the form ( I133P ; such a function is usually called in literature a 
ground state. 

Let us look first at the case a — 1, ip a — ip^ <p a — cp^ for which the equation ( 1141 p yields 
the following representation for G'('0 1 ) from f 1 1 4 1 [) 

G ;te W ) = <^i>M (142) 

Since ip 1 (r) is a monotonic function, we can find its inverse r = r (^ 2 ) , yielding 

(r(s)), O = 0i(oo)<s<^(O). (143) 

Since (r) is smooth, G'(|-0| ) is smooth for < \if)\ < ■0 1 (0). If we do not need G' (s) to 
be smooth, we extend G' (s) for s > ip 1 (0) as a constant, namely 

G\ (a) = Gi (0)) if a > # (0) . (144) 
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G[ (s) 



o 2 

First derivative of such an extension at s = ip x (0) has a discontinuity point. If ip a (r) is a 

° 2 

smooth function of class C n , we always can define an extension of G' (s) for s > ip 1 (0) as a 
bounded function of class C n for all r > and 

G[ (s) = G[ (ft (0)) - 1 if 8 > # (0) + 1. (145) 

Slowly growing functions G' (s) also can be used. 

In the case of arbitrary size parameter a > we define G' a (s) by formula ( 1391) . this 
definition is consistent with (11401) and (11431) . 

Let us take a look at general properties of G' (s) as they follow from defining them 
relations (I143l) - (l39l . In the examples below the G' (s) is not differentiable at s = 0. But if 
•0 (r) decays exponentially or by a power law the nonlinearity g (ip) = G'(\ip\ 2 )ip as it enters 
the field equations is differentiable for all ip including zero, and hence it satisfies the Lipschitz 
condition. For a Gaussian r ip 1 (r) which decays superexponentially G'(\ip\ 2 ) is unbounded at 
zero and g (ip) is not differentiable at zero. Since ip (|x|) > 0, the sign of G[ (\ip\ 2 ) coincides 
with the sign of V 2 ^ (|x|). At the origin x = the function ip 1 (|x|) has its maximum 
and, consequently, G[ (s) < for s close to s — (0). The Laplacian applied to radial 
functions takes the form l x lj- Consequently, if (r) is convex at r = |x| we 

have V 2 ^ x (|x|) > 0. Since r 2 ^p l (r) is integrable, we naturally assume that |x| (|x|) — ► 
as |x| — > oo. Then if the second derivative of rip l (r) has a constant sign near infinity, it 
must be non-negative. For an exponentially decaying ip l (r) the second derivative of rip a (r) 
is positive implying G[ (s) > for s <C 1. In the examples we give below G[ (s) has exactly 
one zero on the half-axis. 

Example 1. Consider a form factor ip l (r) decaying as a power law, namely 

where c pw is the normalization factor, c pw = 3 1 / 2 / (47r) 1 ^ 2 . This function evidently is positive 
and monotonically decreasing. Let us find now G' (s) based on the relations fll43[) . An 
elementary computation of V 2 ?/^ shows that 

^ , . 15s 2 / 5 45s 4 / 5 ^ , > 75s 7 / 5 25s 9 / 5 r a 

° {S) = ^-^ G{s) = ^~^ iOI °- S - C - (147) 

The extension for s > c 2 can be defined as a constant or the same formula (11471) can be 
used for all s > 0. 

If we explicitly introduce size parameter a into the form factor using (I140p . we define 

G'a ( s ) Notice that the variance of the form factor (|x|) decaying as a power law 

ffTtoT) is infinite. 

Example 2. Now we consider an exponentially decaying form factor , 1 of the form 

/ o \i/2 / r°° i o \i/2 \ _1 / 2 

^! (r) = c e e-( r +1 ) , Ce = 4vr / r 2 e" 2 ( r +1 ) dr . (148) 
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Evidently ■0 1 (r) is positive and monotonically decreasing as required. The dependence r (s) 
defined by the relation (I148j) is as follows: 

r = [In 2 (ce/VS) - 1] 1/2 , if Vi < ^ (0) = Cee- 1 . (149) 
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An elementary computation shows that 

-g*L = 2 , +7T l^ + 1 ,-1. (150) 
V>i (r 2 + l)5 (r 2 + l) ( r 2 + i)5 

Combining ( 1149 j) with (j!50p we readily obtain the following function for s < c 2 e -2 

4 4 8 



g; ( S ) 



i - 



(151) 



ln(c 2 / S ) ln 2 (c 2 / S ) ln 3 (c 2 / S ). 

We can extend it for larger s as follows: 

G[ (s) = G[ (c c 2 e- 2 ) = -3 if s > c 2 e" 2 . (152) 

or we can use a smooth extension as in (I145p . The function G[ (s) is not different iable at 
s = 0. At the same time if we set g (0) = the function g (if)) = G[ (if) (r)) if) is continuous and 
g {if)) is continuously differentiable with respect to if) at zero and g (if>) satisfies a Lipschitz 
condition. The variance of the exponential form factor if) 1 (r) is obviously finite. To find 
G' a (s) for arbitrary a we use its representation (1391) . 

Example 3. Now we define Gaussian form factor by the formula 

^{r) = C g e^l\C g = ^ l . (153) 

Such a ground state is called gausson in |Bia] . Elementary computation shows that 

V 2 ^(r) 2 /.a 



r 2 -3 = - In ^ 2 (r)/G 2 ) -3. 



^ (r) 

Hence, we define the nonlinearity by the formula 

G'(|^| 2 )=-ln(|^| 2 /G 2 )-3, (154) 

we call this nonlinearity logarithmic nonlinearity. The nonlinear potential function has the 
form 

G(s) = (- In (s'/G 2 ) - 3) ds' = -s In s + s ^ln - 2^ . (155) 
Dependence on the size parameter a > is given by the formula 



Gl(|^| 2 )=-a- 2 ln(a 3 |^| 2 /G 2 )-3. (156) 

Obviously g (if)) = G[ (\ip\ 2 ) if> is continuous for all if) e C if at zero we set g (0) =0 and is 
differentiable for every t/> 7^ but is not differentiable at if) = and does not satisfy Lipschitz 
condition. 



4 Charges in remote interaction regimes 

The primary focus of this section is to show that if the size parameter a — » the dynamics 
of the centers of localized solutions is approximated by the Newton's law of motion (the 
macroscopic limit a — > assumes that there is a fixed macroscopic scale R ^> a). This 
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is done in the spirit of the well known in quantum mechanics Ehrenfest Theorem, |Schiff| 
Sections 7, 23]. We also Wave-corpuscle solutions defined by (11371) . (11381) provide an example 
of explicit solutions which have such a dynamics. 

As we stated, the Lagrangian £ in (1271) is gauge invariant and every £-th charge has a 
4-current (//, J £ ) defined by 

p = q ^ J = _ _ Im — _ V > 157 

11 \mr tp mrc J 

which satisfies the continuity equations d t p e + V • J e = or explicitly 

dt + V • ( ^ Im^f |^| 2 - ^-A ex |^| 2 ) = 0. (158) 

(Note that J e defined by (I157P agrees with the definition of current (|24p in Maxwell equa- 
tions). Equations fl 1 5 8 [) can be obtained via multiplying fl 1 3 [) by vfj e * and taking imaginary 
part. Integrating the continuity equation we find that ||^|| = const and we impose the 
normalization condition (I133p . The momentum density P e for the Lagrangian £ in dl28[) is 
defined by the formula 

Note that so defined momentum density P e is related with the current J in (I157P by the 
formula 

f 

hi 

P^(t,x) = — J'(t,x). (159) 
We introduce the total individual momenta P £ for £-th charge by 

P^dx, (160) 



and obtain the following equations for the total individual momenta 
dP^ 



dt 

where 



dx, (161) 



v f (t,x) = ip f (t,x) = ij f (^). (162) 

The external EM fields E ex , B cx in (I16ip corresponding to the potentials <£> ex , A ex are deter- 
mined by standard formulas (I20ip . Derivation of (I16ip is rather elementary. Indeed, in the 
simplest case where A ex = we multiply (I130p by V^*, take the real part and integrate the 
result over the entire space using integration by parts. To obtain (11611) in more involved 

general case one can similarly multiply (I130p by V ip * and then integrate the result by parts 
using some vector algebra manipulation. 

Let us show now that if the size parameter a is small compared to the typical scale 
of variation of EM fields, the charge evolution can be described approximately by Newton 
equations with Lorentz forces similar to (jl]). Here we skip technical details. Conditions under 
which this kind of derivation is justified are given in |BabFig3"| . 
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We introduce the £-th charge position r e (t) and velocity v e (t) as the following spatial 
averages: 

v"{t)=v[{t)= I x|^(t,x)| 2 dx, v t (t) = - e [ J^,x)dx, (163) 

where the current density 3 e is defined by (11571) . We show below that a combination of the 
continuity equation (11581) with the momentum evolution equations (11611) imply the following 
remarkable property: the positions r (t) satisfy with a high accuracy Newton's equations of 
motion for the system of N point charges if the size parameter a is small. 

Multiplying continuity equation (I158p by x and integrating we find the following identities 



dt 



[ *d t W\ 2 dx= --. [ xV-j'dx=^/ J £ dx = /(t), (164) 
Jr3 T Jm? Q Jm? 



showing the positions and velocities defined by formulas (11631) are related exactly as in the 
point charge mechanics. Then integrating (1 1 5 9 [) we obtain the following kinematic represen- 
tation for the total momentum 



p i (t) = / t (t, x) dx = mV (t) , (165) 
T Jr3 



which also is exactly the same as for point charges mechanics. Relations (11641) and f II 6 5 j) 
yield 

,dV(t) e d t dP £ 

m = m dt v (t) = IT' (166) 

and we obtain from ( 11611) the following system of equations of motion for N charges: 



,dV (t) 

m 



r > •, , \/ ,.. , - -«v, / |- - • — .n dx, i = 1, N, (167) 



where E e ' (t, x) = —V<p £ ' (t, x), E ex and B ex are defined by (12011) . 

The derivation of the above system is analogous to the well known in quantum mechanics 
Ehrenfest Theorem, [Schifft Sections 7, 23] and |Biaj . Now we give a formal derivation of the 
Newton's law of motion for charge centers. 

Let us suppose that for every E-th charge density \ip \ and the corresponding current 
density J e are localized in a- vicinity of the position r e (t), and that \r e (t) — r e ' (t)\ > 7 > 
with 7 independent on a on time interval [0, T]. Then if a — > we get 

|^| 2 (t, x) -> 5 (x - r e (*)) , v £ (t, x) = J £ /q £ -> (t) 5 (x - i* (*)) , (168) 

where the coefficients before the Dirac delta-functions are determined by the charge normal- 
ization conditions H133[) and relations (I163p . Using potential representations (1 1 3 2 [) we infer 
from (I168P the convergence of the potentials f £ to the corresponding Coulomb's potentials, 
namely 

ip l {t, x) ip' (t, x) = r -^r, V r </ (t, x) 9 lX ~ / as a -> 0. (169) 

|x — r l | |x — r*| 

Hence, if we pass to the limit as a — >■ 0, we can recast the equations of motion (11671) as the 
system 

m^ = f £ + e , (170) 
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where 

f = £ <?X' + 9 £ E ex M + V x B cx (r*) , £ = 1, N, 

and €q — > as a — > 0. Notice that the terms f in equations fjlTOj) coincide with the Lorentz 
forces and we see that the limit equations of motion obtained from (11701) coincide with 
Newton's equations of motion for point charges interacting via the Coulomb forces and with 
the external EM field via corresponding Lorentz forces, namely 

™isr = -£ ^tfz + ^ M + l / x B » M , / = i, - , w. (i7i) 

Note that we essentially use the fact that the nonlinearity G' a , which, according to fl3"5|) . sin- 
gularly depends on a as a — > 0, does not enter the system (j!67p explicitly. For mathematical 



details of the above derivation see BabFig3 



Remark 5 In special case of one particle in an external EM field which depends only on 
time we presented explicit wave-corpuscle solution of U30\) . U31\) given by formula \l3Tty .. 
Note that if the EM fields have a general form of spatial dependence, as they do in the 
case of multiple particles we cannot find write solutions explicitly. Nevertheless, if the EM 
fields do not vary fast in space, which is the case of a system of charges in the regime of 
remote interaction considered in this section, the wave- corpuscle solutions of the form \137\) 
provide approximate solutions to the field equations with high accuracy. The dynamics of 
charge centers r e and their velocities v £ in general case is given by the Newton's equations 
§T7W - Since the spatial extent of if) is small (of order a ) and we can linearize potentials of 
external fields acting on the charge near the center of the wave-corpuscle, the substitution of 
wave-corpuscle into field equations produces relative discrepancy of order 

„ ( a 2 a Ivl 

where R is the typical spatial scale of variation of external fields near the charge trajectory. 
Consequently, we may expect the wave-corpuscle form of localized solutions to be preserved 
for long times. We would like to stress that even for very small a the moving charges are 
not reduced just to points and the oscillatory de Broglie wave factors e ise /x in (Wty remain 
to be significant. More detailed studies of approximate wave-corpuscle solutions in similar 
situation in the framework of the WCM can be found in fBabFigflJ , jBabFig2j. 



5 Hydrogen atom model 

In this section we provide a detailed sketch of our hydrogen atom (HA) model with an 
intention to write a separate detailed paper on this subject. In this model ip e are spinless, 
therefore one cannot expect spin related effects to be modeled here (though it is quite clear 
that extensions of this model to multi-component ip e which can possess spin are possible, this 
is a subject of future research). 

To model the hydrogen atom (HA) we set iV = 2 in the non-relativistic system (11301) . 
(11311) where the indices i take two values £ — 1, for electron, and i = 2, for proton, and the 
charges values q\ = —q = qi- The electric fields in the resting hydrogen atom have to be 
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time-independent, hence \ipe\ 2 in (I13ip must be time-independent too. Therefore we assume 
that only phase factors depend on time and consider the multi-harmonic solutions of this 
system, namely solutions of the form 



V/(t,x) 



V>/ ( X ) 



ip £ (t, x) = <p t (x) 



1,2. 



(172) 



Plugging the expressions (11721) in equations (I130p . f II 3 1 j) we find that the functions ip e (x) 
satisfy the following nonlinear eigenvalue problem 



x 2 



2m 



where, in accordance with (I125p . 



Let us introduce 



— — , 



4^ 



X* 



V 1^ 



q 2 mi ' 



1,2, 



(173) 



(174) 



(175) 



where the quantity a\ turns into the Bohr radius if x equals to the Planck constant H, and 
mi, Q are the electron mass and charge respectively. Using (I175P we rewrite the system (I173p . 
(I174p as the following nonlinear eigenvalue problem 



V 2 $ 2 = -4tt|^ 2 | 2 . 



a 2 , 



a 2 



(176) 
(177) 



V 2 $! = -47r|Vi|, V 2 $ 2 = -4tt|^ 2 | z . (178) 
where ■0 1 and ip 2 are respectively the wave functions for the electron and the proton with 

||^||=1, ||^ 2 || = 1, (179) 



according the charge normalization condition (I133p . The nonlinearities G[, G' 2 are assumed 
to be logarithmic as defined by (I156p where the size parameter a = a 1 is different for electron 
and proton. One can similarly consider other nonlinearities but in this paper we stay with 
the logarithmic ones. 

Let us introduce an energy functional £ ip 2 ) associated with the energy derived from 
the Lagrangian fl 1 28[) by the following formula 



£ (V'u = £1 i>2) + £2 (ipi,ip 2 ) > where 



(180) 



£1 



£2 



q 0,2 



[iV^f + dfl^l 2 )] dx-27rg 2 / (-V 2 )^|^| 2 |^| 2 dx 



[|V^ 2 | 2 + G 2 (|^ 2 | 2 )] dx-27rg 2 / (-V 2 )^|^| 2 |^ 2 | 2 dx 



and (—V 2 ) 1 \ip e \ 2 is defined by the Green's formula (I132p . Observe then that equations 
(11761) . (I177p . (11781) are the Euler equations for critical points of the functional £ under the 
normalization constraint ( I179p . if we set in (I178p 



$! =4tt ((-V 2 ) V 



1 ,.2> 

1 1 ' 



$5 



47T 



(-v 2 )- 1 ^ 
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and the frequencies u}i,u 2 are the corresponding Lagrange multipliers. Importantly, it turns 
out that for the logarithmic as in (11551) nonlinearities Gi, G 2 the frequencies Ui, co 2 and critical 
values of E\,E 2 satisfy Planck-Einstein formula E = ftw (see |BabFig3] for details and also 
[BiaJ where a relation between the Planck-Einstein formula and the logarithmic nonlinearity 
was discovered in a different setting). 

The problem of finding critical values for the energy functional £ defined by (11801) can 
be reduced approximately to a simpler problem for a single wave function in a way similar 
to the Born-Oppenheimer approximation in the quantum mechanics. To that we introduce 
a change of variables 

x=am, i = l,2, (181) 
where ag are defined by (jl75j) . and rescale the fields as follows: 

• t (x)=*M, *,(»)=*£&!>, <=1,2. (182) 
Then the equations (11761) . (11771) turn into the following system 



2 oo 1 ^ 1 + — V^*! + -<j> 2 ( ^yi ) M>i = — G[ (l^l 2 ) * l7 (183) 



a 2 J 


*1 = 


1 

— ( 

2a x 


a 2 \ 

— Y2 

a% J 


*2 = 


1 ( 

2a 2 




h = - 


-4ir 



4^1*1 + + - 

4^2 + 7^V'^ 2 + -0! ( -y 2 ) * 2 = ^G" 2 (|^ 2 | Z ) tt 2 , ( IS I) 

Vj^! = -4tt |^!| 2 , V 2 y2 <P 2 = -47r|* 2 | z , (185) 

where not to complicate notations we use the same letter G in rescaled variables. Recall now 
that the electron/proton mass ratio is small, that is 

&=H*=^JL«i. (186) 
m 2 ai 1800 V ; 

Then we recast (I183p . f !184j) in the following system 



q 2 



+ ^V 2 ^ + Qy) = ^ (l^l 2 ) * l5 (187) 



^2*2 + iv 2 ^ 2 + H x (by) ^2 = \g' 2 (|^ 2 | 2 ) * 2 - (188) 

Notice that equation d 1 8 8 [) for the proton wave function ^ 2 depends on vE^ through the 
potential b<p l (by) which is small since b is small. On the other hand, the electron wave 
functional depends on ^ 2 through the potential j<fi 2 (£y), and if we look for radial solutions 
we can use the formula 



b^ 2 [b r J = r~V J b ~ T ^ n ^ ri ^' 2 dn ' 



(189) 



Restricting ourselves to lower energy levels for £ (■0 1 , i[) 2 ) we can conclude using (11891) that 
\4>2 ii r ) can t> e replaced by the Coulomb potential - with an error of order b 2 (see BabFig3] 
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for details). With that in mind we introduce the following energy functional with the Coulomb 
potential for a single wave function: 



£ Cb i 



<r_ 



±|Wi| 2 + ^(|*i|Vj^| 



dy. 



(190) 



Then with a small and controlled error we substitute the original problem of finding frequen- 
cies Ui based on critical points of £ (V'ijV^) for lower energy levels with a simpler problem 
of finding critical points, lower critical levels and corresponding frequencies U\ for the en- 
ergy functional <?cb (^i) subjected to the constraint H^iH = 1. The corresponding nonlinear 
eigenvalue problem for the electron wave function ^ and dimensionless spectral parameter 



UJ — — 5-Wi IS 



1 



F*l = n G 'l 



2 ) 



(191) 



As we have already mentioned a similar reduction to a single Schrodinger equation with the 
Coulomb potential is made in the quantum mechanics via the Born-Oppenheimer approxi- 
mation. 

Let us exploit the dependence of the nonlinearity G' x = G' la on the small parameter k = — 
which is the the ratio of the electron Bohr radius a\ to the size parameter a. If k is small 
then the nonlinearity G' la (s) = i^G' x (k~ 3 s) is small and plays a role of a small perturbation 
in the eigenvalue problem ( 11911) . A detailed analysis shows that lower energy levels of the 
functional Sch (^i) are arbitrary close to the energy levels of the Schrodinger operator for 
HA provided that k = ^ is sufficiently small. Consequently, based on estimates obtained in 
BabFig3j we can conclude that n-th lower frequency uj\ n for solution of (11761) . (I177j) . (I178[) 



are given by the following approximate formula 



1 q 2 



n 2 2a\ 



a 



ln[^ 



n 



1,2, 



(192) 



The correction term O [b 2 + (^) | In (^) | J in (11921) is small if b given by (11861) and are 

small. Observe that differences of energy levels of the nonlinear eigenvalue problem are very 
close for the same in the Rydberg formula with relative error of order 10 -4 if — is of order 
10 -2 . Hence, if we assume that the size a of a free electron is 100 larger than the Bohr 
radius, then the introduced here hydrogen atom model is a good quantitative agreement 
with the hydrogen spectroscopic data. We think that it is quite reasonable to assume that 
a free electron has much larger size then an electron bound in a hydrogen atom where it is 
naturally contracted by the electric force of the positively charged proton. 

Remark 6 // we use the provisional HA model from [BabFig^ we still obtain discrete 
energy levels, but as — — > the limiting linear eigenvalue problem involves a potential 

2 

+ 90 (y) where in addition to the Coulomb potential there is a term qcf) (y) due the elec- 
tron EM self -interaction. So, if there is EM self -interaction the limiting eigenvalue problem 
as ^- — > does not turn into the same for the linear Schrodinger operator for the HA and 
consequently the energy levels do not converge to the known expressions for HA as — — > 0. 



Now we briefly compare the above non-relativistic treatment of the HA with the treatment 
in the framework of the full relativistic version of our model. We start directly from the 
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relativistic system in (!29|) -( l3T|) and look for time-harmonic solutions with A e = and time- 
independent <fr using substitutions (136]) and (I172p as follows: 

^ (t, x) = e-^+^Hi (x) , <p e (t, x) = Vt (x) , (193) 

2 

where £ = 1, 2, u — = ck m . We arrive then at a system similar to ( I173p . (I174p 



j . - + - ^) 2 ^ + VV ~ G 1 ' (V> V) (—) V = 0. ( 194 ) 

^ V X X J ^ \ X J 

Based on smallness of electron/proton mass ratio we similarly to the non- relativistic case 
arrive to the following eigenvalue problem for electron density which is a relativistic version 
of PU): 

/ 2 \ 2 

(m l c 2 + X ui + ^) ^ 1 + cVV 2 ^ 1 -cVG"(|V 1 | 2 )^ 1 -m2 c V 1 = 0. (195) 

If the ratio k = — is small, the nonlinearity G' (s) = k, 2 G' (k~ 3 s) can be treated as a small 
perturbation, and the linear part of (11951) essentially determines the lower energy levels. Note 
that if we set x — ^ the linear part of (j!95j) coincides with relativistic Schrodinger equation 
(see jSchiff] p. 309 ), which has the form 

2 



(^W 2 + m 2 c>= { E + ]t\) 



u 



with energy level E = mc 2 + x u} i- According to [Schiffj the energy levels of the linear relativis- 
tic Schrodinger equation in a contrast to the non-relativistic hydrogen Schrodinger equation 
have a fine structure, the fine structure energy levels are given by Sommerfeld's formula 
and the relative scale of the fine structure is controlled by a 2 where a is Sommerfeld's fine 
structure constant, a = |— ~ ™. This shows that at atomic scales in our relativistic model 

' nc 137 

relativistic effects are present even in the case of zero velocities if the square of Sommerfeld's 
fine structure constant is not assumed to be negligible. 

Remark 7 In our treatment of charges in Section [7] at macroscopic scales we assume the 
electron size a to be very small and in this section we assume k = ax/a to be very small. 
Since there is huge gap of scales between the macroscopic and atomic scales there is no 
contradiction if we take into account the small value of the Bohr radius a\ ~ 5.3 x 10~ n m 
compared with the scale of variation of EM fields. Note that the error €q of approximation by 
Newtonian trajectory in ^170 ) is of order a 2 /R 2 aaci 1 where R maC r is the scale of spatial 



variation of EM fields which act on a charge. In the treatment of HA in this section we 
assume k 2 = a\/a 2 <^ 1. Taking a ~ lO 2 ^ we arrive at the restriction -R maC r 5.3 x 10~ 9 m 
which is an estimate of the scale of spatial variation of EM fields for which the Newton's 
equations with Lorentz force holds with a good accuracy. 



6 Appendix 

6.1 Classical electrodynamics 

We consider the Maxwell equations for the EM fields and their covariant form following to 
jJacl Section 11.9], jLanLif EMI Sections 23, 30], [GrH Sections 7.4, 11.2], in CGS units 

V • E = Aug, V • B = 0, (196) 
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V x E + -d t B = 0, 

c 



1 47T 

V x B - -dtE = —J. 

c c 



(197) 



To represent Maxwell equations in a manifestly Lorentz invariant form it is common to 
introduce a four- vector potential and a four- vector current density J v : 



(196 



ft, 



g 



c 



<9" 



(9a;,. 



-3t,-V 
c 



and, then, an antisymmetric second-rank tensor, the "field strength tensor, 
so that 






-Ex 


-E 2 


-E 3 


Et 





-B 3 


B 2 


Ei 


B 3 





-Bx 


E 3 


-B 2 


Bx 









Ex 


E 2 


E 3 


Ex 





-B 3 


B 2 


E 2 


B 3 





-B 


E3 


-B 2 


Bx 






and 



E = -Vip - -d t A, B = V x A. 

c 



(199) 



(200) 



(201) 



Then the two inhomogeneous equations and the two homogeneous equations from the four 
Maxwell equations ( I196p take respectively the form 



4:71 

c 



(202) 



d a F M + dpF^ a + d y F aP = 0, a, 0, 7 = 0, 1, 2, 3. (203) 

It follows from the asymmetry of F^", the Maxwell equation ( 12021) and (1198j) - (1199j) that the 
four-vector current J M must satisfy the continuity equation 

8^ = or d t Q + V • J = 0. (204) 

The Maxwell equations (12021) turn into the following equations for the four-vector potential 
A? 

Air 

— J u , (205) 



UA V - d v d ll A tM 



where 



□ = 



— <9 2 — V 2 (d'Alembertian operator). 



The EM field Lagrangian is, |Jacj Section 12.7], |Bar| Section IV. 1] 

4m V^) ~ 



16tt m c M ' 



(206) 



(207) 



where is an external (prescribed) current. Using (I200p . (120 ip and (1 198[) we can recast 
(12071) as 



L em (A") = — (E 2 — B 2 ) 
Sir 



1 

8^ 



V<p + -d t A 

c 



(V x A) 5 



-J A 

c 



pip + - J • A. 

c 



(206 
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In particular, if there are no sources the above Lagrangians turn into 

L em (A") = —F uv F fJ ' u = — (E 2 — B 2 ) = 

v ' 16vr m 8tt v 1 



(209) 



Vp + -d t A - (V x A) 

c 



The canonical energy- momentum (stress-, power- momentum) tensor M " for the EM field 
is as follows, [Jacl (12.104)], jB^l Section III.4.D] 



6 

or, in particular, for i,j = 1, 2, 3 

e 00 = - 
e * = 



F^d v A 



Aix 



1 + 



F&F, 



f7 



167T 



E 2 -B 2 1 <9 A-E 
— + pp J ■ A , 

Sir c 47r 



(210) 



(211) 



diA ■ E 



47T 



, e 



Ekdotp (B x doK), 



47T 



47T 



6 , _ _Ejdj<£_ + (B ^ + g-g _ 1 , . A 

47T 47T 87T C 

whereas the symmetric one a/3 for the EM field is, |Jacl Section 12.10, (12.113)], [Bar, 
Section III. 3] 



~ 4ir V M 4 M 



(212) 



implying the following formulas for the field energy density to, the momentum density g and 
the Maxwell stress tensor r. 



w = @oo = E^ + B 2 ) = @ 0l = Q ,o ^ E x B 



57T 



47T 



6' 



1 

47T 



1 



^ + B^J - -tfy (E 2 + B 2 ) 



(213) 
(214) 



e a/9 



u> eg 



10 
-CE 



-T 



W Cg 
-Cg -Ty 



(215) 



Note that in the special case when the vector potential A vanishes and the scalar potential 
if does not depend on time using the expressions ( 120 1}) we get the following representation 
for the canonical energy density defined by ( 121 ip 



e 



oo 



+ pp for A = and <9o<£> = 0, 

2 



8tt 

whereas 



(216) 



on _ (yvY 
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It is instructive to observe a substantial difference between the above expressions @ 00 , which 
is the Hamiltonian density of the EM field, and the energy density 00 defined by (12131) . 
If there no external currents the with differential conservation laws takes the form 

d a Q a ? = 0, (217) 

and, in particular, the energy conservation law 

1 / dw \ 

= d a Q alB = - I — + V • S , where w is the energy density, and (218) 



c V dt 



S = c 2 g = — E x B is the Poynting vector. 
An 

In the presence of external currents the conservation laws take the form, |Jacl Section 12.10] 

d a Q^ = -f, f = -/^J V) (219) 

and the time and space components of the equations (12191) are the conservation of energy w 
and momentum g which can be recast 

Ht^ s H jE ' (22o) 



dt ^ dxi Tlj 



pEl + I ( J x B) . 

c 



(221) 



The energy conservation law (12201) is often called the Poynting's theorem, |Jac| Section 6.7]. 
The 4- vector f" in the conservation law (12 1 9 [) is known as the Lorentz force density 

f = -F Pv J v = (-3 ■ E, pE + - J x B ] . (222) 
c \c c / 

6.2 Potentials and fields for prescribed currents 

In this section we describe EM fields F^ v arising from prescribed (external) currents J v 
following mostly to |Jac| Section 6.4, 6.5 and 12.11]. Namely, the EM fields F^ u satisfy the 
inhomogeneous Maxwell equation 

d^F"" = ^-J\ F^ = d»A v - d v A», (223) 

which take the following form for the potentials A v 

UA V - d v d lx A^ = -^-r. (224) 

If the potentials satisfy the Lorentz condition, d^A^ = 0, they are then solutions of the 
four-dimensional wave equation, 

OA" = —r (225) 
c 
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The solution of the inhomogeneous wave equation (12251) is accomplished by finding a Green 
function G (x, x') for the equation 



DG (z) = 6 {i) (z) , G (x, x') —G(x- x') 



(226) 



where 5^ (z) = 5 (z ) 5 (z) is a four- dimensional delta function. One can introduce then the 
so-called retarded or causal Green function solving the above equation (J226]), namely 



£(+) ( x _ x >) = 9 (g° I x o) 5 ( Xo - x 'o- R ) R 



(227) 



where 9 (xq) is the Heaviside step function. The name causal or retarded is justified by 
the fact that the source-point time x' is always earlier then the observation-point time xq. 
Similarly one can introduce the advanced Green function 



rH r-r - ^ - 9 t~ ( x ° ~ x ' ^ 5 ^° - x 'o + R ) p 

Lr [X X ) - ^ R ' 



x — X 



These Green functions can be written in the following covariant form 

1 



G {+) (x - x') 



2tt 



9 (X — Xq) 5 



x — x 



(228) 



(229) 



G { - ] {x-x') = —9(x' -x )5 (x-xf 



where 



x — x) 2 = (xq — x' Y — |x — X' 

= 2R [ tf ( x o ~x' -R) + S {x -x' + R)} 



I \2 



J\2 



(X — X 



The more explicit form of the Green functions G^ in terms of time space variable is 



where 



R 



T = t-t', 



or 



The solution to the wave equation 



:S t 



(230) 



(231) 



(232) 



(233) 



can be written in terms of the Green functions 



47T 



A v (x) = A\ xl {x) + — {x - x') J v {x') dx' 



or 



A" (x) = A" out (x) + 



An 



G^ (x - x') J v (x') dx' 



(234) 



(235) 



where A\ n (x) and A^ ut (x) are solutions to the homogeneous wave equation. In ( 12341) the 
retarded Green function is used. In the limit x — > -co, the integral over the sources 
vanishes, assuming the sources are localized in space and time, because of the retarded 
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nature of the Green function, and A\ n (x) can be interpreted as " incident" or " incoming" 
potential, specified at x$ — > — oo. Similarly, in (I235j) with the advanced Green function, the 
homogeneous solution A^ ut (x) is the asymptotic "outgoing" potential, specified at xq — > +00. 
The radiation fields are defined as the difference between the "outgoing" and "incoming" fields, 
and their 4-vector potential is, [DirGE] 

47T f 

A^ d (x) = A» out (x) - Al (x) = —jG(x- x') r (x') dx', where (236) 

G(x- x') = G {+) {x - x') - (x - x') . 

More explicit form of the potential A u (x) solving the inhomogeneous wave equation (I225P 
based on the retarded Green function and with A" Q (x) = is 

(t> x) = f [pgijOLt dx / 5 A (t, x) = - / [J dx', (237) 

J R cj R 

where 

R = x-x', |x-x'|, (238) 

and the symbol [-] ret means that the quantity in the square brackets is to be evaluated at the 
retarded time 

R Ix — x'l 
t' = t rct = t =*-■! l -. (239) 

c c 

The corresponding to the potentials I237f) EM fields can be represented by Jefimenko formulas, 
|Jef[ Section 15.7], |Jacj Section 6.5] 



E (t, x) = / Lrv ; ^ Jret dx' + / ^ /Jret dx'- (240) 



n^j(t', x ')] rct 
y c 2 i? d ' 



B(( ,x)^/{Ii(^ + M|^}xR dx ', ,24!) 

where R = y§i[- ^ n essentially equivalent form of the Jefimenko equation (I240p for the electric 
field E is due Panofsky and Phillips, [PanPhit Section 14.3] 

E(f ,x) = / ^'g^ ft dx' + (242) 
, v [J (f, x')]«, • A) A + ([J (f , x')L„ x A) x R 

> t .J(i',x')] rel xR) xR 

It was pointed out by McDonald in [McDoJ that the combination of equations ( 124 ip and 
( 12421) has a certain advantage since it "manifestly displays the mutually transverse character 
of the radiation fields (those that vary as 1/R)". Since the radiation fields E ra d and B ra d 
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decay asl/R for large R we can extract them from the expressions (12421) . (12411) for the entire 
EM fields obtaining 



E rad (t, x) 



([^J(t',x')] rct xR) X R 
^ =r= '- dx' 



c 2 R 



B rad (t,x) = J [^J(t',x')] ret x ^dx' 
6.3 Point charge and the Lienard-Wiechert Potential 



(243) 
(244) 



If the particle is a point charge q whose position and velocity in an inertial frame are respec- 
tively r (t) and v (t) = d t r (t) the corresponding charge and current densities in that frame 
are, [Jacl Section 12.11, (12.138)] 



p (t, x) = qS (x - r (t)) and J (t, x) = qw (t) 5 (x - r (*)) . 



(245) 



Using the formulas (123 7p for the charge density and the current as in the equations (12451) we 
obtain the Lienard-Wiechert Potential, [Jac[ Section 14.1] 



(p (t, x) 



1-/3R 



A(t,x) 



ret 



qf3 



1-/3R 



J ret 



where 



R = x-r(t r ), R=\x-r{t t )\, R = (3 = ^ 

R c 



and the retarded time t T = t T (x, t) is defined implicitly by the following equation 



(246) 



(247) 



(248) 



Then with the help of the Jefimenko formulas (I240p , (124 ip applied for the point charge density 
and the current (I245p one can derive the Heaviside-Feynman formulas (first discovered by 
Heaviside (1902) and rediscovered by Feynman (1950)), [Heavlj . |Heav2[ Subsection 510], 
|Fey[ Vol. I, Section 28; Vol II, Section 21], |Jac[ Section 6.5], [Jan], |Monaj : 
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(249) 
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R 



x d t 



R 



+ 



ret 



E (t, x) x 



R 








ret x d 2 


R 


c 






R 








ret 





- ret 
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where 



x= 1 



v R 



(251) 



In view of the implicit relations ( 1248ft between the retarded time t T and the time-space 
variables (x, t) it is important to keep in mind that there is evidently a difference between 
dt [(OU and [d t (-)] ret . 

The Heaviside-Feynman formulas (I249p . (I250P imply the following formulas for the radi- 
ation fields of the moving point charge 



E ra d (t, x) 



^d 2 

2 * 



R 



ret 



B ra d X) 



C 2 [fl] 



-a 



ret 



v x R 



(252) 



Another important representation of the EM fields of an arbitrary moving charge is their 
decomposition into the velocity and acceleration fields, |Jacl[ Section 14.1] 



E (t, x) = E v (t, x) + E a (t, x) , where for (3 



E v (t, x) = q 



R-(3) (1-/3 2 ) 



3 R 2 



J ret 



E a (t,x) = E rad (t,x) = ~ 



R 



x 



R-/3) x /3 



- ret 



(253) 



B(i,x) = RxE (t, x) = B v (t, x) + B a (t, x) 
B v (i, x) = R x E v (t, x) , 
B a (t, x) = B rad (t, x) = R x E a (t, x) . 



(254) 



The velocity fields are essentially static fields falling off as R~ 2 , whereas the acceleration 
fields are typical radiation fields, both E a and B a being transverse to the radius vector R and 
varying as R~ l . For low velocities the formulas ( 12531) . ( 1254}) turn into the following simpler 
asymptotic expressions 



Erad (t, x) 



B r ad (t, X) 



q 



Rx R x v 



R 



(255) 



ret 



v x R 
IRI 



(l + 0(\(3\)), 



ret 



6.4 Almost periodic functions and their time-averages 

We provide here very basic information on real and complex valued almost periodic (a. p.) 
functions following |Cort Section 3.1]. In fact, we consider for simplicity sake a class A\ (R, C) 
of a. p. function of the form 

oo 

/ (*) = J2 ^ e " iaJs< ' -oo < t < oo, (256) 

s=l 
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where the exponents u s are real valued numbers and amplitudes f s are complex valued 
numbers and such that 

oo 

i/i = X)i/.i< 00 - ( 257 ) 

8=1 

We refer to the set of exponents Af = {oo s } as the Fourier spectrum of the function / and 
to the numbers f s as its Fourier coefficients. The class Ai (R, C) is a Banach algebra, that is 
(i) it is linear space and (ii) for any / and g in A\ (R, C) the product fg is in A\ (R, C) as 
well and \ fg\ < \ f\ \g\. The derivatives of a.p. f (t) satisfy 



dU(t) = J2(-^sYfse-'^ t , (258) 



s=l 



provided 



J2\"s\ r \fs\ <oo. (259) 



8=1 



Every a.p. function / in Ai (R, C) is assigned the time average (mean) value (/) defined by 

(/>= ^2?// (t)di - <260) 
The mean value has the fundamental property 

</<«> •*> = {£' olT^e • < M > 

In some applications it is convenient to view a real a.p. function / (t) as the real part of 
a complex valued function f (t) written in the following form 

f (t) = Yl /- e ~^ (262) 

where f u = (f (t) e iujt ) = (f cos (cot)} + i (/ sin (ut)) . 

The set Af in f)262p is at most countable set of non-negative frequencies u > 0, and we refer 
to it as the frequency spectrum of /. Then the corresponding real valued a.p. function has 
the following representation 

/ (t) = Re {f (t)} = \ (f.e' iut + ft**) • (263) 

Evidently, any real a.p. function can be represented in the form ( 12621) . (12631) . Observe that 
for any a.p. complex valued function f and g of the form f)262p we have 

| V Re{U*J ifA / nA,^0, 
<Re{f(t)}Re{g(t)}> = { weA/ nA 9 • (264) 

ifA/nA 9 = 

In particular, if / and g are real valued the relation (I264j) imply 

| Y\ Re{U*J ifA / nA 9 ^0, 
(f(t)g(t)} = { uieAftiAg ■ (265) 

ifA/nA 9 = 
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The identity (I264p readily follows from formula f l262[) . It reads that if the frequency spectra 
and A+ have no common frequencies the time-average (Re {/} Re {g}) is identically zero. 

A particular case of the formula (I264p when the frequency spectra of / and g consist of a 

the same single frequency u is customary used in electrodynamics for time harmonic fields, 

[PanPhil Section 11.2], [Strati Section 2.20]. 

For any a.p. function / with the frequency spectrum Af one can introduce the smallest 

additive group in the set of all real numbers that contains all frequencies u from Af. Such a 

smallest group is called the module of / and is denoted by mod (/), |Cor| Section 4.6]. It is 

easy to see that mod (/) consists of all real numbers of the form 

s 

rrijCUj where rrij are integers, Uj E Af and s is a natural number. 
6.5 Vector Identities 

Here is the fist set of commonly used vector identities 

a ■ (b x c) = b • (c x a) = c • (a x b) , (266) 

ax (b x c) = (a-c)b- (a • b) c, (267) 
(a x b) ■ (c x d) = (a ■ c) (b • d) - (a ■ d) (b ■ c) . (268) 
Using the above identities we readily obtain for any vectors a, b and any unit vector u 

u x (a x u) = a — (a • u) u, |u| = 1, (269) 



[u x (a x u)] x (u x b) = [(a x u) x u] x (b x u) = (270) 
= [a • b — (a ■ u) (u • b)] u = 
= [(a x u) ■ (b x u)] u, |u| = 1. 



If 2 x 

— / [a-b-(a-x)(x-b)] da = -a • b, x = — . (271) 
4tt J| x |=i 3 |x| 

V • (a x b) = b • (V x a) - a • (V x b) . (272) 
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